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RISPOSTA DINAMICA DELLE STRUTTURE 



1.1 L 'Oscillatore semplice 

Una struttura che la cui una massa pud schematizzarsi come concentrata in un punto e delta 
"oscillatore semplice". Segue da tale schema che lo stato di deformazione strutturale e noto 
quando e noto il moto della massa, unica fonte delle forze di inerzia che agiscono sulla struttura. 
Esempi di oscillatore semplice sono mostrati in Figura 1. 
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Figura 1 esempi di strutture schematizzabili come un oscillatore semplice e schema strutturale 

Con riferimento alia figura si osserva che le masse del ponte e del serbatoio sono essenzialmente 
concentrate in sommitd. Nel ponte si osserva che le travi tampone dell 'impalcato sono ciascuna 
appoggiate da un lato su un carrello e dall'altro su una cerniera, pertanto il ponte in senso 
longitudinale pud essere schematizzato con tanti oscillatori semplici (se si assume nullo I 'attrito dei 
carrelli) costituiti ciascuno da una pila priva di massa (massa trascurabile) e dalla massa in 
sommitd, somma di quella della stampella e della trave tampone collegata alia stampella con una 
cerniera. 
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Le due strutture indicate in Figura 1 possono quindi essere schematizzate come una massa 
sopportata da strutture elastiche di sostegno prive di massa. Nella figura si e indicato lo schema a 
portale con trove infinitamente rigida. La rigidezza alia traslazione orizzontale e K (K/2 per 
ciascun piedritto). 

E' possibile scrivere I'equazione di equilibrio alia traslazione della massa. Detta : 

x,x,x } Q spostamento relativo della massa rispetto alia base e le sue derivate prima e seconda 

rispetto al tempo (velocitd ed accelerazione relativa) ed ? la accelerazione della base rispetto ad 
un sistema fisso, cioe I 'accelerazione del terreno, la forza di inerzia Fi e data dal prodotto della 

massa m per I 'accelerazione assoluta x + ^ : 
F t =m-(x + y) 

La forza elastica di richiamo nella posizione di quiete: x=0 vale 
Fe e = k-x 

Vi e poi una forza dissipativa, di tipo quindi non conservativo, che e comodo aifini degli sviluppi 
per la soluzione elastica in forma chiusa mettere nella forma di dissipazione viscosa: 
F d =d-x 

Infine pud agire sulla massa una generica forza F. 



Per I 'equilibrio della massa si deve avere che ad ogni istante: 
F l+ F d +F e =F 

sostituendo si ottiene: 



m -('x + y) + d ■ x + k ■ x = F 



(1 



Naturalmente nella espressione (1) la accelerazione del terreno y e la storia della forza 
esterna F sono note, oltre ad essere note m, k, d, caratteristiche meccaniche della struttura, 
pertanto si pud po rtare a secondo membro i termini noti: 



m-x + d- x + k- x = F- m- y 



(2 



i termini m- x;m- y hanno le dimensioni di una forza ma non hanno significato fisico, la loro 

somma e la forza di inerzia. Tuttavia e consuetudine chiamare il secondo "forza di trascinamento " 
della struttura. 
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1.1.1 Soluzione dell'equazione di equilibrio 

Risolvere I'equazione differ enziale ( 2 vuol dire trovare la funzione x(t) che soddisfa I'eguaglianza. 
La ( 2 e una equazione differenziale a coefficienti costanti. La soluzione si trova come somma 
dell' integr ale generale, che e la soluzione dell'equazione ( 2 con termine noto nullo, che 
rappresenta I'equazione di equilibrio in assenza di moto del terreno e di forze esterne ed e detta 
equazione omogenea associata, V integrale particolare che soddisfa il caso con il termine noto 
assegnato, nel caso in questione F -m-y . 
Si divida la (2 per la massa m, si ottiene I 'equazione: 



.. d 


. k 


F .. 


(3 


x + — 


x + — 


■x = y 


m 


m 


m 





si vede quindi che le tre grandezze m, d, k non intervengono nel I 'equilibrio in maniera 
indipendente tra loro, ma legate da un rapporto; si pone allora 



k 2 d d 


(4 


— = co ; — = 2vco ^>v = 


m m 2am 





v e detto rapporto di smorzamento, (essendo il rapporto tra il coefficiente di smorzamento d ed il 
prodotto 2a>m detto smorzamento critico) e vale per le strutture reali tra lo 0.5 ed il 5%, 
quest 'ultimo valor e essendo quello tipico degli edifici. 

Come si vedrd nel seguito, una struttura avente uno smorzamento pari al critico, e quindi v=l, 
caso che non si verifica per alcuna struttura dell'ingegneria civile, allontanata dalla sua posizione 
di equilibrio, vi ritorna asintoticamente senza compiere alcuna oscillazione. 



1.1.1.1 SOLUZIONE PER IL CASO DI TERMINE NOTO UGUALE A ZERO 
(oscillazioni libere) 



X + 2VG) -X + CO ■ X = 0 



(5 



per v = 0 L 'integrale generale, che rappresenta la soluzione del problema, essendo il termine 
noto nullo, e del tipo: 



x = Acoscot + Bsincot 



(6 
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in cui le costanti A e B devono essere determinate ponendo le condizioni ai limiti, in questo caso si 
ha per t =0 x = x 0 ;x = x 0 (condizioni iniziali) : 



x () = Acos(a>0) + Bsin(a>0) =^> x 0 = A 

X 

x 0 = -coAsin(coO) + <yficos((yO);=> x 0 = coBcos(coO) => 5 = — 

CO 



sostituendo i valori trovati nell'integrale generate si ha: 



x = x 0 cos cot H — sin cot 

CO 



L'espressione della risposta mostra che la storia degli spostamenti di una struttura priva di 
smorzamento e periodica, essendo somma di due funzioni periodiche che assumono lo stesso valor e 
dopo per valori dell 'argomento cot che differiscono di un multiplo di 2n. Pertanto ragionando in 
termini di tempo la risposta assume gli stessi valori quando: 



cot = 2n => t = — = 2n 

CO 




II valore T viene pertanto detto periodo proprio della struttura, ed e una grandezza che dipende 
dalle caratteristiche meccaniche della stessa. Esso e il tempo che intercorre tra due oscillazioni 
consecutive di massima ampiezza di una struttura nelle sue oscillazioni libere. Si noti che esso 
aumenta con la massa e diminuisce con la rigidezza. I valori del periodo delle strutture civili 
variano in genere tra qualche decimo di secondo ed al piu qualche secondo. Ad esempio per gli 
edifici a telaio vale la regola empirica: T=0.1N ove N e il numero dei piani, col che si capisce che 
in Italia la gran parte degli edifici ha periodi compresi tra 0.1 e I secondo (1-10 piani). Vengono in 
genere dette "rigide" le strutture con periodo inferiore a 0.2-0.3 secondi, "medie" tra 0.3 e 0.6- 
0.7 secondi flessibili al di sopra di 0.8-1 secondo. I terremoti in generale impegnano 
particolarmente gli edifici con periodo compreso tra 0.1 e 0.6-0.8 secondi, cioe proprio il campo di 
periodi propri dell 'edilizia corrente italiana. Sono meno sensibili i grattacieli, che hanno periodi 
spesso superiori ai 3 secondi ed i grandi ponti, si pensi ad esempio che il Golden Gate Bridge di S 
Francisco ha il periodo proprio in senso trasversale di 18.2 secondi ed in senso verticale di 10.2 
secondi. 



7 



Prof. Camillo Nuti Dispense di Costruzioni in Zona Sismica 2007 



Per 0< v <l, (le strutture reali hanno sempre v >0, in genere v non supera il 5-10%) la soluzione 
e: 



x = e vmt [A cos co x t + B sin co x t] . 






(9 


x = e [x 0 cos co x t + — sin co x t\ ; 




CO 




co x = coyjl-v 2 





Vespressione tra parentesi e analoga a quella trovata per v = 0 con la sola differ enza la frequenza 
co e leggermente minore a causa dello smorzamento. Questa diminuzione della frequenza e 
trascurabile per i valori di smorzamento tipici delle strutture civili: si tenga presente che se 
v = 0.20 coi=0.98co. 

Si hanno pertanto oscillazioni di frequenza e quindi di periodo pressoche uguali a quelli pre- 
cedenti, ma di ampiezze decrescenti per effetto del termine esponenziale negativo che contiene il 
rapporto di smorzamento v. 

Analogamente a quanto fatto in precedenza si determinano le costanti A e B imponendo le 
condizioni iniziali e si ottiene Vespressione ( 10: 



vwtr X n +VCOX n . 

x = e [x 0 cos co J + — -sincoj] 

co x 

;co x = co^jl-v 2 


(10 


X 

x _ e ~v<* CQS (0 f + _0_ 
CO 

y 


(11 
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J.52U 
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Figura 2 Funzioni di risposta x(t,v), per diversi valori dello smorzamento. La funzione x(.) e 
valutata con Vespressione approssimata ( 11, la funzione xl(.) con Vespressione esatta ( 10. I 
valori adottati sono m=25.35 t, k=1000 kN/m. Spostamento iniziale Xo =3 m, velocita iniziale =40 
m/sec. Si noti Vinfluenza di v. 



II prodotto vcaxo e trascurabile rispetto a x 0 , cosi come a>i differisce poco da co, pertanto dal punto 

di vista pratico e possibile utilizzare Vespressione della risposta non smorzata moltiplicata per il 
termine smorzante esponenziale, secondo Vespressione (11, come peraltro dimostrato dalla Figura 
2. 

Da questa espressione sembrerebbe che una struttura piii rigida (cogrande) a parita di 
smorzamento v, smorzi le oscillazioni piu rapidamente; ma in realta questo non avviene spesso 
poiche non e esattamente valido che la viscositd sia direttamente proporzionale alia velocita, come 
e stato assunto nello scrivere la equazione del moto. 

v=l (smorzamento critico) manca il termine oscillatorio e si ha: 



x = e vm '[x 0 (l + (Ot)] 



(12 



e un caso che non presenta nessun interesse pratico. 

Esaminiamo il caso avente condizioni iniziali : spostamento nullo e velocita assegnata. 
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x 0 = 0 => x 0 = A = 0 

X 

x 0 = coBcos(coO) => B = — 

CO 



La legge del moto diviene quindi: 



x = e L — sinco x t\ 

CO, I-, 2 

1 ;co l = coyll-v 


(13 


x = e [ — sin cot \ 

CO 


(14 
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Figura 3 Funzioni di risposta x(t, v), nel caso di spostamento iniziale Xo =0, velocita iniziale =40 
m/sec. La funzione x(.) e valutata con Vespressione approssimata ( 13 la funzione xl(.) con 
Vespressione esatta ( 14. 1 valori adottati sono m=25.35 1, k=1000 kN/m. 
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II problema fisico e lo studio di una vibrazione provocata da una forza del tipo F(t) applicata per 
un tempo estremamente limitato. Si definisce impulso la seguente espressione: 



/ = lim {F(t)dt 

A£->0 J 




x(t) 






\ / \ / t 
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del moto Riscrivendo I'equazione 



m'x + k ■ x = F(t) 



(16 



si possono integrare vari membri da to to+As calcolarne il limite per Asche tende a 0: 



tn+S t 0 +E 



lim \mxdt+ \k-xdt = \F(t)dt 

\e^>0 J J J 

? o l o *o 

lim f m — dt+ \k-xdt = \F(t)dt 
l<r->0 J At J J 

t t i 

'o 'o 'o 

t 0 +E t 0 +E t a +£ 

lim \mdx+ \k-xdt = \F(t)dt 



'0 '0 




ss 
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Ottenendo, per II termine noto, Vespressione dl una forza impulsi va. 

II primo termine, In cui si e espressa 1 ' accelerazione come derivata della velocita rispetto al 
tempo, rappresenta, dopo la semplificazione, I 'integrate di un differenziale di una funzione, quindi 
la funzione stessa. II secondo termine invece si pud risolvere solo esplicitando la funzione x(t), ma 
essendo lo intervallo infinitesimo possiamo approssimarla con una funzione linear e x(t)= at 
Si ottiene allora: 



m[x]' 0+£ + k ■ a[—]' a+s = I 

o 2 ° 



il primo termine e la differenza tra la velocita finale e quella iniziale moltiplicata per la massa, ma 
essendo la velocita, la velocita iniziale e tuttavia nulla, pertanto essa rappresenta la velocita che ha 
la massa su cui ha agito Vimpulso, e quindi e la velocita del moto che ne deriva. II secondo termine 
al tendere dell 'intervallo a zero, risulta essere un infinitesimo di ordine superior e e quindi 
trascurabile; si ha quindi: 





■ T ■ 1 

mx Q = I => x 0 = — 
m 


(19 




Sostituendo quindi la legge del moto derivante da una forza impulsiva assume 
espressione: 




x = ——e va " sin cot 
m co 


(20 




chiamata "funzione di risposta a impulso e Indicata con il simbolo h(t). 




Se I 'impulso anziche all'istante t=0 agisce all'istante t=r, allora lo spostamento x 
dovuto all 'impulso dipende dal tempo intercorso At=t-r, pertanto si ha 




X {t) = /(r) 1 e - wo{ '- T) sinco(t r) 
m co 


(21 
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La generica storia diforze F(t), pud essere pensata come successione di impulsi F( r)dr 
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Figura 4 Esempio diforza continua 
ed effetto all'istante 4.5 di tre 
impulsi per t=l,2,3 e loro somma (in 
basso a destra) 



f(t) o 




La risposta al tempo t diviene quindi la sovrapposizione degli effetti di tutti gli impulsi che 
precedono il tempo t, cioe: 



m a> 



sinoo(t - r)dr 




Assumendo come F( r) proprio il termine noto della equazione di equilibrio (2, e sostituendo nella ( 
22 si ottiene: 



. . \-ma(r) 1 - va> ( t -r\ . , , 

x(t) = — — e ( } sinoo(t - r)dt 

* m co 

x(t) = \a{T)—e- vco{, - T) sinco(t-T)dT 

{ CO 



(23 
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La espressione ( 23 e detta risposta ad un accelerogramma nel dominio del tempo o anche 
"integrale di Duhamel". 

Volendo eseguire un calcolo approssimato, ma piu efficiente dal punto di vista del tempo di calcolo, 
si pud osservare che e sufficiente tener conto degli impulsi che precedono Vistante t di calcolo di un 
tempo multiplo del periodo di non oltre 10 volte, ad esempio per una struttura con T=0.5, e 
sufficiente considerare un tempo precedente di non oltre 5 secondi. 
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1.2 Lo spettro di risposta 



Osservando la espressione ( 23 si nota come, dato un accelerogramma, la risposta dipenda 
esclusivamente dai due parametri, ved co. Pertanto due strutture aventi gli stessi valori di ved co 
hanno la stessa storia della risposta, ed in particolare la stessa risposta massima. In generale nelle 
costruzioni civili non si e interessati all'intera storia della risposta ma ai valori massimi delle 
sollecitazioni e degli spostamenti. Una volta noto lo spostamento massimo Xmax, la sollecitazione 
massima vale: 




E ' possibile cosi dimensionare la struttura se si vuole che resti in campo elastico. In alternativa, 
noto Xmax, e possibile dimensionare la struttura perche sopporti lo spostamento Xmax. 
Dal punto di vista del progettista e pertanto molto utile dato un accelerogramma disporre di un 
diagramma dove in ascisse sono riportati i periodi propri T ed in ordinata il valore di Xmax. Si 
possono costruire curve per diversi valori del rapporto di smorzamento v. 

Si osserva che ricavata lo spostamento massimo, e la relativa forza, e possibile chiedersi quale e la 
accelerazione che produrrebbe una uguale forza di inerzia: 



F m »v = k ■ X m „ = m- Sa 

max max 



(25 



Invertendo la ( 23 si ottiene: 





(26 


Sa = --X m , ix ^Sa = co 2 -X mi 


m 





Se tutta I'energia di deformazione si trasformasse in energia cinetica, cosa che avverrebbe in 
assenza di smorzamento, si avrebbe: 



E = \-k-X 2 mm . = ±.m-Sv 2 mm ^ Sv max = co ■ X n 




Dalla (23 e (23 si ottiene il legame tra Sa ed Sv: 
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Sa =co- Sv 



(28 
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Ec8 Amax=1(m/sec A 2) 
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Figura 5 Spettri di risposta in spostamento, pseudovelocita e pseudoaccelerazione, per 
diversi valore dello sporzamento. In basso spettro di risposta in scala trilogaritmica. 
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Figura 6 spettro di risposta in velocita in scala logaritmica. 
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Sv ed Sa sono dette Pseudovelocitd e Pseudoaccelerazione deU'oscillatore, non essendo 
esattamente la velocitd o I'accelerazione assoluta dello stesso. Si osserva in particolare che la 
accelerazione assoluta deve tendere al valore della accelerazione del terreno per T che tende a 
zero, cioe per strutture infinitamente rigide che traslano quindi con il terreno, inoltre la velocitd 
relativa deve tendere alia velocitd del terreno quando la struttura e infinitamente flissibile, quando 
cioe T tende aU'infinito, poiche la massa in tal caso rimane ferma rispetto al sistema fisso mentre 
il terreno si muove. Vi sono pertanto deifferenze tra pseudovelocitd e velocitd relativa della 
struttura in corrispondenza dedi periodi lunghi, mentre la accelerazione assoluta differisce dalla 
pseudoaccelerazione in prossimitd di T=0. 

Si osserva infine che, come detto,nel caso di struttura molto flessibile, la massa tende a rimanere 
ferma, pertanto lo spostamento relativo e pari a quello assoluto del tereno Uq, pertanto lo spettro 
di risposta in spostamento deve tendere a U per T che tende aU'infinito. 
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1.2. 1 Spettri di risposta di accelerogrammi 



Spettri di risposta su roccia 



(a) 



"6 



7 eventi sismici 
(14 accelerogrammi) 




San Rocco5 NS 
San Rocco5 EW 

Tarc4 NS 
Tarc4 EW 

Robic3 NS 
Robic3 EW 
Gebze2 NS 
Gebze2 EW 
Hercegnovi2 NS 
Hercegnovi2 EW 
Robic4 NS 
Robic4 EW 
San Rocco3 NS 
San Rocco3 EW 
Media 

EC8 - Sott. A 



Periodo, T (s) 



d 

N 
I 



10 eventi sismici 
(20 accelerogrammi) 




1 2 
Periodo, T (s) 



Gebzel NS 
Gebzel EW 

Hercegnovil NS 
Hercegnovil EW 

Ulcjnil NS 

Ulcjnil EW 

Bagnl NS 

Bagnl EW 

Sturnol NS 

Sturnol EW 
S.G. La Molara NS 
S. G. la Molara EW 
Tolm2 NS 
Tolm2 EW 
Hercegnovi5 NS 
Hercegnovi5 EW 
San Rocco4 NS 
San Rocco4 EW 
Robicl NS 
Robicl EW 

Media 

EC8 - Sott. A 



Figura 7 spettro di risposta per accelerogrammi reali a confronto con il loro spettro medio 
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e lo spettro normativo su roccia per terremoti intensi (Tipo 1) 

Gli spettri di risposta sopra illustrati si riferiscono ad accelerogrammi registrati su roccia 
selezionati (vedi Rey et al. 2002) 1 dallo European Strong Motion Database, con particolare 
riferimento a quelli relativi ai terremoti del Friuli, Montenegro e Izmit. Gli accelerogrammi sono 
stati normalizzati dividento le ordinate per il picco di accelerazione, pertanto Vordinata dello 
spettro di risposta al periodo proprio T=0 vale 1. 

Una operazione che spesso viene effettuata per ricavare spettri di progetto e quella di fare la 
statistica di spettri di risposta rappresentativi, con qualche criterio, della situazione di interesse. 
Una analogia di minimo e quella delle caratteristiche di intensitd epicentrale, e di situazione 
geologica e geotecnica del sito ove si vuol valutare lo spettro. 

Nelle figure sono mostrati anche gli spettri ottenuti mediando quelli relativi ai diversi 
accelerogrammi. Si ottengono cost forme piu regolari di quelle relative ai singoli accelerogrammi. 
Si noti come a intensitd epicentrale maggiori (in magnitudo), riportati nella figura inferior e, 
corrispondano accelerogrammi con spettri di risposta con la zona di massima amplificazione piu 
estesa. Questa situazione e rilevabile in modo pressoche sistematico. Pertanto e da attendersi che 
gli spettri su roccia relativi a zone di pericolositd sismica inferior e, cioe quelle ove Vazione sismica 
attesa e minore, dovrebbero avere una estensione della zona di massima amplificazione minore. 



Rey J. , Faccioli E., Bommer J. (2002). Derivation of design soil coefficients (S) and response spectra shapes for eurocode 8 using the European 
Strong Motion Database. Journal of Seismology, 547-555 
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1.2.2 Spettri di risposta ricavati a partire dalle caratteristiche di sismicita regionale 

Dato V epicentro di un terremoto di nota intensitd sono state messe a punto per via statistica 
funzioni, dette leggi di attenuazione, che a partire dalle caratteristiche di intensitd epicentrale, 
danno, in base alia distanza del sito ove si vuole valutare lo spettro di risposta, il valore delle 
ordinate spettrali medie. Con opportune operazioni di media si possono cosi calcolare, a partire 
dalle caratteristiche di sismicita locale, spettri di risposta di sito, in genere su roccia, sempre su 
terreno a superficie orizzontale, che tengono conto delle caratteristiche di sismicita del territorio 
circostante. 

A titolo di esempio si riportano i valori ricavati per alcune cittd Italiane [Nuti, Rasulo & Vanzi, 
2005J 2 . 



2 Nuti,C, Rasulo, A., Vanzi, I., (2005) "Proposta per la Valutazione della Classificazione Sismica del Territorio Italiano " Rapporto Tecnico Nl/2005 
Dipartimento di Strutture Universita degli Studi RomaTre. 
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CATANIA 



MESSINA 



NAPOLI 



Messina, Zona 1°, Catania e Napoli, Zona 2°, Confronto tra spettri di risposta T=475 anni 
ottenuti a partire dalla recente classificazione INGV, ottenuti con la classificazione GNDT 
2001 dal SSN, Rivalutati da Sabetta Pugliese da [Nuti, Rasulo& Vanzi 2005] 
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0,45 



0,40 



0,35 



0,30 



0,25 



0,20 



0,15 



0,10 



0,05 



0,00 



NutiRcSuloVanzi 




ServizioSismipo 



RW Romeo 



0,00 



0,50 



1,00 



1,50 



2,00 



2,50 



3,00 



3,50 



4,0C 



Confronto tra spettri di risposta con periodo medio di ritorno 475 anni ottenuti a partire dalla 
recente classificazione INGV, ottenuti con la classificazione GNDT 2001 dal SSN, Rivalutati da 
Sabetta Pugliese, da [Nuti, Rasulo& Vanzi 2005] 



La tecnica piu diffusa per la valutazione statistica a partire dai dati di sismicitd del territorio e 
quella mess a a punto da Cornell nel 1969. Per un approfondimento si rimanda a testi piit 
specializzati. 

Le forme spettrali cosi ottenute dipendono oltre che dalle leggi di attenuazione dalla definizione 
delle caratteristiche della sismicitd regionale. Si noti ad esempio che gli spettri ottenuti con la 
recente classificazione dell 'INGV 2004, sono piu bassi di quelli ottenuti a partire dalla precedente 
classificazione del GNDT 2001. 

Si noti come le forme spettrali cosi ottenute sono anch 'esse piuttosto regolari, e presentano una 
certa analogia con quelle ottenute con operazioni di media di spettri di accelerogrammi naturali. 
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La cosa e d 'altronde da attendersi in quanto le leggi di attenuazione sono ricavate per via statistica 
dalle registrazione degli accelerogrammi naturali. 

Se invece degli spettri con periodo medio di ritorno 475 anni si valutano quelli con periodo 72 
anni, corrispondenti a probabilitd di superamento del 50% in 5° anni, si possono ottenere forme 
spettrali quali quelle della figura seguente: 



0,14 



0,08 



0,06 
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Catania 
Firenze 

Genova 

Messina 
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Venezia 

Verona 


I 

I 


I 


























































































































































































































































































i 




I 

































































































































































0,00 0,50 1,00 1,50 2,00 2,50 3,00 3,50 4,00 4,50 



Spettri medi con periodo medio di ritorno 72 anni. 
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1.2.3 Spettri di risposta risposta di norma 

Le normative nazionali ed internazionali definiscono in genere le forme spettrali in relazione alle 
caratteristiche delle formazioni piu superficiali del terreno. Ad esempio I'Eurocodice 8 fa 
riferimento alle caratteristiche degli strati sino a 30 metri di profonditd. Questa profonditd non ha 
un particolare valore scientifico, piuttosto rappresenta un valore sino al quale si spingono i 
sondaggi che si effettuano per la valutazione delle caratteristiche geotecniche dei terreni di 
fondazione quando si deve costruire un edificio. II fondamento scientifico si basa sulla valutazione 
che spesso a 30 metri di profonditd gli strati di terreno divengono molto compatti per cui non ha 
interesse pratico indagare le caratteristiche a profonditd piu grandi. 

Qualora le caratteristiche del terreno siano di qualitd scadente a profonditd maggiori, e pertanto 
opportuno fare studi specifici. 

E' tuttavia ormai riconosciuto che le forme spettrali attese in un sito specifico dipendono non solo 
dalle caratteristiche degli strati superficiali, ma anche dalle caratteristiche alia sorgente sismica e 
dai terreni attraversati. Pertanto si ritiene che sia quantomeno opportuno che le forme spettrali 
vengano definite a livello regionale attraverso studi specifici che tengano conto di tutti i possibili 
parametri, non ultimo la morfologia locale. 



L ' Eurocodice 8 da la seguente espressione per lo spettro elastico di progetto: 



o<r<r fl :5 e (r)=a 8 -5 



l+— (ti-2,5-1) 



r B <7 <7 t :S e (r)=a g -5-Ti-2 ) 5 



7 C < 7 < 7 n : S e (T) = a g S r\ 2£ 



T D <T< 4s: 5 e (r) = a St] -2,5 



T C T D 
T 2 



n =^10 /(5+4)>0,55 



smorzamento spettrale 
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Settri elastici per terremoti intensi: Tipo I 
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Ground type 


S 


Ms) 


7'r (s) 


'/'d(s) 


A 


1,0 


0.15 


0.4 


2,0 


15 


1.2 


0.15 


0,5 


2,0 


c 


1.15 


0.20 


0.6 


2,0 


L) 


1,35 


0.20 


0,8 


2.0 


L 


1.4 


0. 1 5 


0.5 


2.0 



Valo ri dei parametri nell 'EC8, terremoti intensi: Tipo I 
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Settri elastici per terremoti deboli: Tipo II 
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Ground type 


s 


7b (s) 


7c (s) 


7n(s) 


A 


1.0 


0.05 


0.25 


1,2 


13 


1,35 


0.05 


0.25 


1.2 


C 


1,5 


0. 1 0 


0.25 


1.2 


L) 


1.8 


0.10 


0.30 


1.2 


1. 


1.6 


0.05 


0.25 


1.2 



Valo ri dei parametri nell'EC8, terremoti deboli: Tipo II 
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La proposta di annesso tecnico nazionale dell'EC8, di gennaio 2005, rporta i seguenti valori dei 
parametri 



Zona sismica 


Categoria 


S 


TB 


TC 


TD 




suolo 












A 


1,00 


0,10 


0,40 


4,50 


1-2 


B 


1,15 


0,15 


0,60 


5,00 




C-D-E 


1,30 


0,20 


0,80 


6,00 




A 


1,00 


0,05 


0,35 


1,50 


3, 4 


B 


1,20 


0,10 


0,45 


1,50 




C-D-E 


1,35 


0,15 


0,60 


2,00 



Valori dei parametri dello spettro di risposta elastico della componente verticale di accelerazione 
per lo SLU 



Zona sismica 


Categoria 


S 


TB 


TC 


TD 




suolo 












A 


1,00 


0,10 


0,40 


4,50 


1-2 


B 


1,15 


0,15 


0,60 


5,00 




C-D-E 


1,30 


0,20 


0,80 


6,00 




A 


1,00 


0,05 


0,35 


1,50 


3, 4 


B 


1,20 


0,10 


0,45 


1,50 




C-D-E 


1,35 


0,15 


0,60 


2,00 



Valori dei parametri dello spettro di risposta elastico delle componenti orizzontali di accelerazione 
per lo SLD 



Lo spettro di risposta relativo alia componente verticale e generalmente concentrato su frequenze 
piu elevate di quello per le componenti orizzontali. 
L Eurocodice 8 da le seguenti espressioni: 

o<r<r B :s ve (/-)=« Ji+^--(n -3.0-1) 



T B <T<T C : S ve (r)=a vg Ti.3,0 
T C ZTZT D : 5 ve (r)=a vg ri-3,0 



T D <T <4s: S W (T)= a -r)-XO 



/ parametri sono dati nella seguente tabella: 
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Spectrum 


a va /a a 


r B (s) 


T c (s) 


r D (s) 


Type 1 


0.90 


0.05 


0.15 


uo 


Type 2 


0,45 


11 ( J S 


0.15 


1.0 



Nell 'annesso tecnico nazionale si danno i gli stessi valori. 
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1.2.4 Confronto tra spettri normativi ed alcune proposte in lettratura 



Confronti EC8 / Annex Italiano / OPCM 3274/PitiIakis 
Terremoti Tipo 1 
Sottosuolo Tipo A Sottosuolo Tipo B 



A-EC8 
A - Annex Naz. 
A - OPCM 
A - Pitilakis 




5 

B 

i 
CJ 
C- 

t/3 

u 

B 
O 

1 



B -EC8 
B - Annex Naz. 
B - OPCM 
B1,B2- Pitilakis 




Periodo, T (s) 



Periodo, T (s) 
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Sottosuolo Tipo C 



Sottosuolo Tipo D 



C-EC8 

C1,C2- Annex Naz. 
B,C,E - OPCM 
C2,C3 - Pitilakis 




1 2 3 

Periodo, T (s) 



E-EC8 

C6 - Annex Naz. 
B,C,E - OPCM 
E - Pitilakis 




73 

i= 



D-EC8 

C3,C4 - Annex Naz. 
D - OPCM 
D1,D2 - Pitilakis 




1 2 3 

Periodo, T (s) 



Sottosuolo Tipo E 



1 2 3 

Periodo, T (s) 
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Confronti EC8 / Annex Italiano / OPCM 3274/Pitilakis 
Terremoti Tipo 2 

Sottosuolo Tipo A Sottosuolo Tipo B 




01234 01234 
Periodo, T (s) Periodo, T (s) 



Sottosuolo Tipo C Sottosuolo Tipo D 




0 12 3 4 

Periodo, T (s) 



Prof. Camillo Nuti Dispense di Costruzioni in Zona Sismica 2007 



1.2.5 Testo della Ordinanza della Protezione Civile 3274 relativo alia definizione della 
Azione Sismica 

Si riporta il testo dell' allegato tecnico alia Ordinanza di Protezione Civile 3274. Si tratta di una 
formulazione simile a quella dell'EC8 e relativi Annessi tecnici nazionali. Per questi ultimi si veda 
anche il sito Web: www, coordinatore. it . 

3. AZIONE SISMICA 

i.\ Calegorie di suolo di fondazinnt' 

Ai fini della definizione della azione sismica di progetto si definiscono le seguenti categorie di profilo stratigrafico del suolo di 
fondazione (le profondita si riferiscono al piano di posa delle fondazioni. i valori da utilizzare per Vs. NSPT e Cu sono valori 
medi): 

A - Fonnazioni litoidi o snoli omogenei molto rigidi caratterizzati da valori di V s30 superiori a 800 m/s. comprendenti 
eventuali strati di alterazione superficiale di spessore massimo pari a 5 m. 

B - Deposit! di sabbie o ghiaie molto addensate o argille molto consistenti, con spessori di diverse decine di nietri. 
caratterizzati da un graduale miglioramento delle proprieta meccaniche con la profondita e da valori di V sl0 compresi tra 
360 m/s e 800 m/s (ovvero resistenza penetrometrica N spT > 50, o coesione non drenata c>250 kPa). 

C - Depositi di sabbie e ghiaie mediamente addensate. o di argille di media eonsistenza, con spessori variabili da diverse 
decine fino a centinaia di metal, caratterizzati da valori di V s io compresi tra 180 e 360 m/s (15 < N spr < 50. 70 <c u <250 
kPa). 

D - Depositi di terreni gramtlari da sciolti a poeo addensati oppure coesivi da poco a mediamente consistenti , caratterizzati 
da valori di Vyj < 180 m/s (N spr < 15. e u <70 kPa). 

E - Profili di terreno costituiti da strati superficiali allitvionali. con valori di Vyo simili a quelli dei tipi C o D e spessore 
compreso tra 5 e 20 m. giacenti su di un substrato di materiale piti rigido con V s i 0 > 800 m/s. 

In aggiunta a queste categorie. per le quali nel punto 3.2 vengono definite le azioni sismiche da considerare nella 
progettazione. se ne definiscono altre due. per le quali sono richiesti studi speciali per la definizione dell'azione sismica da 
considerare: 

51 - Depositi costituiti da. o che includono. uno strato spesso almeno 10 m di argille/limi di bassa eonsistenza. con elevato 
indice di plasticita (PI > 40) e contenuto di acqua. caratterizzati da valori di V S30 < 100 m/s (10 < c u < 20kPa) 

52 - Depositi di terreni soggetti a liquefazione. di argille sensitive, o qualsiasi altra eategoria di terreno non classifieabile nei 
tipi precedenti 

Nelle definizioni precedenti V lJ0 e la velocita media di propagazione entro 30 m di profondita delle onde di taglio e viene 
calcolata con la seguente espressione: 




(3.1) 



dove hi e V, indicano lo spessore (in m) e la velocita delle onde di taglio (per deformazioni di taglio y < 10" b ) dello strato i- 
esimo. per un totale di N strati presenti nei 30 m superiori. 

II sito verra classificato sulla base del valore di V sM . se disponibile. altrimenti sulla base del valore di Nspt- 



36 



Prof. Camillo Nuti Dispense di Costruzioni in Zona Sismica 2007 



3.2 C'iilcnlo cli-ll'iiziniu- MMiika 

.12.1 Zone sismiche 

Ai fini dell'applicazione di queste norme. il territorio nazionale viene suddiviso in zone sismiche. ciascuna contrassegnata da 
un diverso valore del parametro a g = accelerazione orizzontale massima su suolo di categoria A (definito ai punto 3.1). con 
probability di superamento del 10% in 50 anni. I valori di a g , espressi come frazione dell'accelerazione di gravita g. da adottare 
:n ciascuna delle zone sismiclie del lenilorio nazionale bono, salvo piti accurate deleiininazioni. clie pussono poilaie a 
differenze comunque non superiori al 20% dell'accelerazione per le zone 1 e 2 e non superiori a 0.05g nelle altre zone: 



Zona 


Valore di a s 


1 


0.35g 


; 


0.25g 




0.1 5g 


4 


0,05a 
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3.2.2 Descrizione dell'azione sismica 



II modello di riferimento per la descrizione del moto sismico in un panto della superficie del suolo e costituito dallo spettro di 
risposta elastico di cui al punto 3.2.3. Quulora siano eseguite determinazioni piii accurate del moto sismico atteso. 
corrispondenti alle probability di superamento definite in 2.1 e 2.2. e consentito utilizzare spettri specifici per il sito pure he" le 
ordinate di tali spettri non risultino in nessun punto del campo di periodi di interesse inferiori all'80% delle ordinate dello 
spettro elastico standard (p. 3.2.3) applieabile in relazione alia categoria di suolo (p. 3. 1 ). 

Per applicazioni particolari. il moto del suolo pu6 essere descritto mediante accelerogrammi. secondo quanto indicato al punto 
3.2.7. 

II moto orizzontale e considerato composto da due component] ortogonali indipendenti. caratterizzate dallo stesso spettro di 
risposta. 

In mancanza di documentata informazione specifica. la componente verticale del moto sismico si considera rappresentata da 
uno spettro di risposta elastico diverse da quello delle componenti orizzontali. come specificato in 3.2.3. 

3.2.3 Spettro di risposta elastico 

Lo spettro di risposta elastico e costituito da una forma spettrale (spettro normalizzato), considerata indipendente dal livello di 
sismicifa. moltiplicata per il valore della accelerazione massima (agS) del terreno che caratterizza il sito. 

Lo spettro di risposta elastico della componente orizzontale e definito dalle espressioni seguenti: 

0<T< T B S e (T) = a g S I 1+— (n-2,5-l)| 

S.(T) = a.S r, 2.5 (3.2) 



T C <T<T D S.(T) = a, S n 2.5|^ 

T C T D 



S,(T) = a. S n. 2,5 



T 



nelle quali: 

S fattore che tiene conto del profilo stratigrafico del suolo di fondazione (vedi punto 3. 1 K 

Per sbuttuie con fattore di importanza yf> 1 . di cui al paragrafo 4.7. erette sopra o in vicinanza di pendii con inclinazione > 15° 
e dislivello superiore a circa 30 metri. fazione sismica deU'equazione (3.2) dovra essere incrementata moltiplicandola per un 
coefficiente di amplificazione topografiea S T . In assenza di stndi specifici si potranno utilizzare per S T i seguenti valori: 

a) S T = 1.2 per siti in prossimita del ciglio superiore di pendii scoscesi isolati; 

b) S T = 1.4 per siti prossimi alia sommita di profili topografici aventi larghezza in cresta molto inferiore alia larghezza alia 
base e pendenza media > 30°; 

c) Sj= 1.2 per siti del tipo b) ma con pendenza media inferiore. 
II prodotto S*Srpu6 essere assunto non superiore a 1.6. 

f\ fattore che tiene conto di un coefficiente di smorzamento viscoso equivalents % diverso da 5 01=1 per 5=5). 
essendo c espresso in percentuale : 

il = v i0/( 5 + ;) > 0,55 (3.3) 
T periododi vibrazione delfoscillatore semplice; 

T B , T c , T D periodi che separano i diversi rami dello spettro. dipendenti dal profilo stratigrafico del suolo di 
fondazione. 

I valori di T B . T c . T D e S da assumere. salvo piti accurate determinazioni. per le componenti orizzontali del moto e per le 
categorie di suolo di fondazione definite al punto 3.1 . sono riportati nella Tabella 3. 1. 
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Tabella 3.1 - Valcri del parametri nelleespresslonl (3.2) dello spettro dl risposta elastico delle componemi orizzontali 



Categor.a suolo 


S 


T B 


T c 


T„ 


A 


1.0 


0.15 


0.40 


2.0 


B. C. E 


1,25 


0,15 


0.50 


2.0 


D 


1.35 


0.20 


0.80 


2.0 



Lo spettro di risposta elastico delta componente verticale e definito dalle espressioni seguenti: 
0<T<T B S ve (T) 

T B <T<T C S W (T) 
T c < T < T D S ve (T) 

T D <T S ve (T) 

con i valori dei parametri che definiscono la forma spettrale riportati in tabella 3.2. 

Tabella 3.2 - Valori dei parametri dello spettro di risposta elastico della componente verticale 



Categor.a suolo 


S 


T B 


T c 


T D 


A. h C. I). 1. 


. i 1 


i i.i i : : 


0.15 


1 J' 



1 ii spettro di risposta elastico dello spostanienlo potril ottenersi jer trasforma7ione diietta dello spettro di risposta elastico delle 
atcelerazioni. usando la segue nte espressione: 

S De (T) = S e (T) ^| (3.5) 

Gli spettri sapra defimti potrantu essere applicati per periodi di vibrazione che non eccedono 4,0 s. Pet periodi superiori lo 
spettro dovra essere definito da appositi studi. 

Nei casi in cui non si possa valutare adeguaiamente rappartenenza del protllo stratigrafico del suolo di fondazione ac una delle 
categorie di cui al punto 3. 1. ed escludendo comunque i profili ci tipo SI e S2, si adcttera in generale la categoria Do. in caso 
di incertezzt di attribuzione tra due categorie. la condizione piO cautelativa. 

3.2.4 Spostamento e vdocitd del wrreno 

I valori dello spostamento e della velocita orizzontali massimi del suolo (rf s ) e (v 8 ) sono dati dalle seguenti espressioni: 
d£= 0.025 S Tr Tn a E 

(3.6) 

v 8 =0,16S T c a s 

3.2.5 Spettri di progetto per lo stato limite ultimo 

Ai ftni del progetto. le capacita dissipative delle strutture possono essere messe in conto attraverso tin fattore riduttivo delle 
fcrze elastiche. denominato fattore di struttura q. L'azione sismica di progetto S d (T) e in tal caso data dallo spettro di risposta 
elastico di cui al punto 3.2.3. con le ordinate ridotte utilizzando il fattore q. I valori numerici del fattore q vengono defimti in 
Rmzione dei matenali e delle tipologie stmtturali. come indicato successivamente nelle presentt norme. 
Lo spettro di progetto per le componenti oiTzzontali e definito dalle seguenti espressioni: 



0.9a S- 1+ (n 3.0-1) 



= 0.9a c S V| 3.0 



(3.4) 



=0.9a g S 



=0.9a g S 
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0<T<T B 
T„ <T<T C 
T C <T<T D 



r t (25 

S d (T) = a g -S- 1 + — ■ 

S d (T) = a 0 S — 



S d (T) = a e S — (-^ 

q 1 t 



(3.7) 



S (1 (T) = a It S 



23 i TcJd. 
q I T 2 



in cui Tb, Tc. Td sono definiti in tab. 3.1. Si assumera comunque Sj(T) > 0.2a g . 

A meno di adeguate analisi giustificative, lo speltro di progelto della componente verticale dell'azione sismica e dato dalle 
seguenti espressioni. assume ndo q = 1.5 per qualunque tipologia strutturale e di materUle: 



O _ T T u 



T„ <T< T r 



S vd (T) = 0.9a S 



S vd (T) = n.9 a, S 



S vd (T) = 0.9-a 8 -S 



S rf (T> = 0.9 a £ S 



T 3.0 

t r I q 



3.0 



(3.8) 



30| Jc_ 

q t 



3,0 I T r T D | 



q I T 



in cui Tfj, Tc- Tp sono definiti in tab. 3.2. Tutti i simboli mantengono significato e valore numerico definiti nel caso dello 
spettro di risposta elastico 



3.2. f> Spettro di progetto per lo state limit? di danno 



Se non si esegue una puntuale valutazione dell'azione sismica corrispondente alia probability di superamento di cui al punto 
2.2. lo spettro di progetto da adottare per la limitazione dei danni puo essere ottenuto riducendo lo spettro elastico di cui al 
punto 3.2.3 secondo un fattore pari a 2,5. 
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1.3 Oscillatore Semplice: Soluzione nel caso in cui il termine noto e una forza 
cosinusoidale 

L'espressione d"equilibrio della massa del sistema, in assenza di smorzamento, e la seguente: 



mx + d ■ x + k ■ x 


= F Q cos fit 








(29 


x + 2vco ■ x + co 2 


■x = —cosilt 






m 





Ove D rappresenta la frequenza della eccitazione esterna. In questo caso la soluzione 
dell 'equazione differ enziale e la somma dell 'integrale generale (uguale a quello trovato nel caso di 
oscillazioni libere) e dell 'integrale particolare che varia nel tempo con la stessa frequenza 
della eccitazione esterna indipendentemente dalla struttura: 



x = e vmt [A cos a> x t + B sin cof\ + [A x cos Qt + B sin Q?, ] 



(30 



Nella espressione ( 23 il primo termine della addizione e I 'integrale generale, esso all 'aumentare 
di t tende a zero, la legge del moto a regime, quindi, sard rappresentata dal integrale particolare: 



x(t) = [A { cos Q,t + 5, sin £lt x ] 



(31 



Basterd allora calcolarsi I valori delle costanti Ai e Bj. Si deriva quindi due volte I'integrale 
particplare ricavando le espressioni di della velocitd e accelerazioni relative in funzione delle 
costanti A i e Bj incognite. Queste espressioni, sostituite ai vi valori nella equazione differ enziale, ci 
permettono di calcolare lare i valori diAjeBj che la soddisfano. 



1 m <V-Q 2 ) 2 + 4v 2 <yQ 



s, - — 



2vaQ. 



m (co 2 -Q 2 ) 2 + Av l o£l 



(32 



Sostituendo queste espressioni nell'integrale particolare si ha la soluzione: 
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x{t) = ^ 



m ^(a, 2 _ Q 2 ) 2 + 4v 2 coQ. 
F 

x (t) = -5- AT (v, co, Q) • cos(Q? - (/>) 
m 



cos(Q? - <j>) 



M(v,co,0) = 



1 



^{oo 1 -O 2 ) 2 + 4v 2 ooQ 



= ^A 2 + B 2 



. ,B.. 2va£l 
<p = artg(—) = 



4 (<y 2 -Q 2 ) 




Ove la fonzione M( v,co,Q) e detta funzione di trasferimento, (/> e I'angolo di sfasamento della 
risposta. 



Si vede, quindi, che I'ampiezza della risposta (Z direttamente proporzionale all "intensita della 
eccitazione ma dipende anche, tramite la funzione di trasferimento dalla frequenza dell 
'eccitazione, da quella della struttura, dallo smorzamento. 
Esplicitando la funzione M(.) si ha: 



M(v,oo,Q) = - 

Jd-— ) 2 + 4v 2 — 
V co co 






CO 


(34 


M(v,co,n) = -L 1 = J_ . M{y ,p) 

<° \J(1-J3 2 ) 2 + 4v 2 J3 2 or 




»(y,P)= i — — 

j(l-j3 2 ) 2 +4v 2 /3 2 





sostituendo questa nuova funzione nell ' espressione della risposta ( 33 si ottiene: 



x(t) = — ~- //(v, (3) ■ cos(Q^ -(/>) = 
m co 

F m 

x(t) = -5- — • //(v, J3) ■ cos(Qt -(/>) = 
m k 



(35 
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F 

x(t) = — ■ ju(v, /?) ■ cos(Clt -(/)) = x u ■ ju(v, P) ■ cos(Clt - (/)) 
k 




7.5 



H(0.01,P) 



\i (0.05, P) 




lO.124, 0 I 1 =i 

0 0.5 1 1.5 2 2.5 3 

0 p A 



Figura 7 Andamento del coefficiente di amplificazione /4v,P) per diversi valori dello 
smorzamento v. 



Dalla terza delle ( 35 si pud dare una interpretazione delta risposta dinamica. Essa e data 
essenzialmente dalla stessa espressione della risposta statica, rapporto tra forza e rigidezza k, 
sfasata rispetto alia applicazione della forza di un tempo t=<f)/Q, ed amplificata di un fattore pari a 
ju( v,P) detto per questo "funzione di amplificazione dinamica". 



La funzione di amplificazione dinamica dipende dallo smorzamento e dal rapporto tra la frequenza 
eccitatrice e quella della struttura. 

Se P=0 si ha il caso di forza statica (0=0) o di struttura infinitamente rigida (co=co). II valore di 
/u(v,p=0) vale 1. Se /3=<x>, si ha il caso struttura molto flessibile con periodo proprio molto lungo 
(co=0) o di eccitazione ad alta frequenza (Q=co) (a>=<x>). il valore di ju( v,p=co)=0. Se (3=1, si ha il 
(co=Q) cioe di frequenza propria uguale alia frequenza di eccitazione, il valore di 
/u(v,f3=l)=l/(2v), pertanto se lo smorzamento e piccolo il coefficiente di amplificazione diviene 
molto grande. La situazione e detta di risonanza, ed al limite, dopo molte oscillazioni la struttura 
pud amplificare moltissimo I 'effetto dell 'azione. 
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1.3.1 Valutazione del transitorio 



Se si vuole valutare il transitorio si deve risolvere I'equazione differ enziale ( 23. L'integrale 
particolare dell 'equazione del moto e quello illustrato, si ricava d'altronde sostituendo il secondo 
termine della soluzione generale, cioe la soluzione particolare nella equazione di equilibrio, 
ricavando le espressioni di A e B sopra illustrate: 



x = e va * [A cos co x t + B sin mf\ + [^4 1 cosQ? + 5 1 sinQ? 1 ] 


(36 




i _F 0 co 2 -a 2 

1 m (a 2 -Q 2 ) 2 + 4v 2 an 
B _ F 0 2vo£l 
1 m (co 2 -Q 2 ) 2 + Av 2 ojQ, 


(37 



Derivando la espressione ( 23 si ottiene I ' espressione della velocitd, si possono cosi imporre le 
condizioni iniziali e da queste ricavare le espressini di A e di B in funzione di Au B\ e dello 
spostamento e velocitd iniziali: 



x Q = [Acos coft + B sin cofi] + [A l cosQ0 + B { sinQOJ 

x Q = -vcoe~ y<m [A cos co x t + BsincoJ] + e~ vwt [-Aco x sin&y + Bco x cos co x t] + [-AflsinQt + .SQsinQ^] 


(38 




Xq — A ~\~ A^ 

x 0 = -vcoiAcoscofi] + Bco x cos^O + ^QcosQjO] = -vcoA + Bco x + Bfi 


(39 




A — JCq A^ 

x t) vco t B.Q. x n vco , JS B,Q. 
B = — + — A + — — = — + — (x 0 -A l ) + - } — 
co l co x co l co l co x co x 


(40 
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-0.129,0.2 



0 
0 



0.5 



1.5 

t 



2.5 



Figura 8 Forzante cosinusoidale. In alto: risposta nel caso diperiodo proprio T=0.127sec, 
periodo delta forzante 0=2 sec, coeff di smorzamento v=0.05, ampiezza delta forzante 
F=100, massa 0.4 (la rigidezza e quindi pari a 1000).In basso: risposta nel caso di periodo 
proprio T=0.28 sec e massa 2.0, gli altri parametri sono immutati. 



x: risposta totale, xl: integrale generale, x2: risposta senza contributo del transitorio 

Le espressioni ( 40 sono diverse da 0 anche nel caso in cui velocitd e spostamento iniziale sono 
nulli, pertanto il transitorio da sempre luogo ad una alterazione della risposta, che tende ad 
annullarsi con le oscillazioni. Ad esempio, in Figura 8 si nota che la risposta totale x(t), e diversa 
inizialmente dalla sola risposta a regime x2(t), la differenza e data dal contributo del transitorio 
xl(t). 

Nella figura e ben evidente il fatto che I 'applicazione della forzan te gid diversa da 0 da luogo ad un 
considerevole differenza nella fase iniziale tra forzante e risposta, differenza che tende a ridursi 
dopo poche oscillazioni. L 'effetto del transitorio e piii modesto se viene applicata una forzante del 
tipo Fsin(fk) che parte anch 'essa dal valore nullo al tempo t=0. 
Nel caso di forzante sinusoidale le costanti A] e Bj divengono: 
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m (co 2 -Q 2 ) 2 + 4v 2 a£l 



F 0 2va£l 

A, — — 



m (co 2 -CI 2 ) 2 +4v 2 6)Q. 




Sono pertanto uguali alle espressioni della forzante cosinusoidale ma sono scambiate ed una, la 
A i con segno invertito, mentre le A e B dell'integrale generate si ottengono imponendo che la 
soluzione completa soddisfi le condizioni iniziali, hanno pertanto le medesime espressioni ( 40, ove 
A j e Bi hanno le espressioni appena tr ovate. 
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0.2 

,.0.103 j 

0.1 

x(t) 

xl(t) 0 
x2(t) 

0.1 

.-0.102,0.2 

( 
























































0 0.5 1 1.5 2 2.5 
) t 


3 
3 


Figura 9 Forzante sinusoidale: risposta net caso di periodo proprio T=0.28sec, periodo 
delta forzante &=2 sec, coeff di smorzamento v=0.05, ampiezza delta forzante F=100, 
massa 2 (la rigidezza e quindi pari a 1000) (tranne la forzante e it secondo caso delta 
Figura 8) 

x: risposta totale, xl: integrale generate, x2: risposta senza contributo del transitorio 
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1.4 Oscillatore semplice: Risposta al moto del terreno 

Se il moto del terreno e del tipo A(t)=ao sin cot, I'equazione di equilibrio e la stessa ( 29 nella quale 
a termine noto e il termine -m ao sin cot, anziche la forzante F sin cot. La risposta e quindi quella 
appena vista ove al termine F si sostituisce -m ao- II rapporto tra la accelerazione assoluta del 
sistema e la accelerazione alia base, e detto tramittanza del sistema: 



(x(t) + a(t)) r 



( 



l + (2vw/Q) 2 



(1-O/Q) 2 ) 2 +2v(co/n) 2 



)' 



(42 



Questa espressione e la medesima che si ottiene nel caso di forzante sinusoidale tra forza totale di 
risposta e forza agente: 

(dco-x + k-x) mm 



La stessa relazione vale tra spostamento assoluto, somma di quello relativo e quello del terreno, e 
spostamento del solo terreno. La trasmittanza ha pressoche lo stesso valore del coefficiente di 
amplificazione, salvo per valori molto elevati dello smorzamento. 



10 



10 



7.5 



Tr(0.01,fi) 



Tr(0.05,n) 
Tr(0.10,Q) 
Tr(0.20,n) 



2.5 



.0.042, o 



0 
0 















\ 








If 


I 

i 
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60 



80 
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Figura 10. Andamento delta trasmittanza in una struttura con frequenza propria w=10 per 
diversi valori del coefficiente di smorzamento 



Si noti che la risp osta diviene quindi: 



x(t) 



ma n 



yjico 2 -Q 2 ) 2 +4v 2 coQ 



m 



cos(Q£ - (/)) 



x(t) = a 0 M(v,co,Q.)-cos(Q.t-tf>) 



(43 
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M(v,co,n) = ] —= =L = ^A 2 +B 2 

■J (co 2 -Q 2 ) 2 + 4v 2 coQ 



. S lx IvcoCl 

<p = artg(—) 



A' {co 2 -a 2 ) 



La funzione M() e detta funzione di trasferimento, se si vuole la risposta in accelerazione la 
funzione va derivata due volte rispetto al tempo, ottenendo quindi la funzione di trasferimento in 
accelerazione relativa: CfM(). 



M(v,co,Q) = - 1 

Jd-— ) 2 + 4v 2 — 
V co co 




CO 


(44 


M(v,co,n) = ~ 1 = \ ■ ju(v,/3) 
(o ^](\-j3 2 ) 2 + 4v 2 j3 2 o) 




yj(l-j3 2 ) 2 +4v 2 /3 2 
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1.5 Oscillatore semplice: Valutazione delta Risposta in Forma Numerica. 

1.5.1 Calcolo della risposta dinamica al passo: Metodo di integrazione al passo 
dell'equazione del moto 

In genere, disponendo di un calcolatore, e conveniente integrare numericamente I'equazione di 
equilibrio dinamico: 



mx{t) + dxit) + kx{t) = -ma(t) (1 



piuttosto che non risolvere I'integrale di Duhamel o la soluzione relativa alia sovrapposizione in 
frequenza. Si tratta di discretizzare il tempo in intervalli di dimensione costante At, sufficientemente 
piccoli e scrivere I'equazione dinamica non in forma differenziale ma in forma finita. La funzione 
spostamento e le sue derivate sono calcolate in corrispondenza dei tempi ti, in cui e stato 
discretizzato il fenomeno. 

I metodi numerici rappresentano I 'unica soluzione per il calcolo della risposta se la struttura entra 
in campo non lineare. Peraltro stime della risposta non lineare possono essere fatte con metodi 
approssimati basati su approssimazioni lineari della risposta. 

Si illustrano nel seguito due metodi di integrazione numerica: il primo delle differenze centrali, nel 
quale I 'incognita di spostamento all'istante i+1, e espressa in funzione dello spostamento e delle 
sue derivate (velocita ed accelerazione) agli istanti precedenti: i ed i-1, e detto per questo un 
metodo esplicito; il secondo, chiamato metodo di Newmark, esprime lo spostamento all'istante i+1 
in funzione degli spostamenti agli istanti i ed i-1, ma anche della accelerazione all'istante i+1, lo 
stesso in cui viene calcolato lo spostamento, e che e anch 'esso una incognita. II metodo, almeno in 
linea di principio, richiede di iterare a partire da un valore di tentativo della accelerazione 
all'istante i+1, ed e per questo detto implicito. 



1.5.2 Metodi espliciti: le differenze centrali 

La velocita all'istante i pud essere espresse come rapporto incrementale tra quelle agli istanti i+1 
ed i-1: 

■ _ x i+\ ~ x i-\ 



L 'accelerazione si pud ottenere dal rapporto tra I'incremento della velocita tra i+1 ed i , ed i ed i- 
1: 
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x. = ■ 



At 



At 



X, = 



X i+\ 2x ; + X t _ v 



(Atf 



(3 



(4 



Sostituendo la (3 e la (2 nella equazione di equilibrio dinamico (1 si ottiene: 



X;,-\ X: 



m 



A ' M 1V ' ' V;| • ' • A.v, = ma(t : ) 



(Atf 



2At 



(5 



Raccogliendo i termini con lo stesso indice e portando a destra dell 'uguale i termini che dipendono 
da i ed i-1 si ottiene: 



m d 

+ 



(A?) 2 2 At 



x M =-ma{t i ) 



m d 




j 2m 


(Atf 2At_ 







(6 



La (6 e una equazione ricorrente, all'istante i+1 sono note tutte le grandezze relative agli istanti 

precedenti, pertanto e possibile ricavare lo spostamento x relativo ad 

La (6 pud ess ere vista come la classica Kx=F, in cui la F e il secondo membro. 

Si noti inoltre come nella (6 compaia il termine Kx,, che rappresenta la forza elastica all'istante i. 
Qualora la struttura abbia comportamento non lineare I 'espressione relativa alia forza di richiamo 
e f(x), pertanto la (1 diviene : 

mx(t) + dx(t) + / (x(t)) = -ma(t) (7 

e la (5 diviene: 



m d 
[(Atf ' 2At_ 


x M =-ma(t i )- 


m d 
_(Atf 2At_ 


x M - 


2m 

^ — (a/) 2 _ 




(8 



Si deve osservare come anche la forza di richiamo non lineare sia nota in quanto e proprio quella 
relativa al passo precedente. Nessuna differenza, nemmeno di complessitd di calcolo, vi e in questo 
caso nel passaggio da una analisi lineare ad una non lineare. 



Si deve infine osservare che per il calcolo del primo passo t=At=i+l, non e possibile definire lo 
spostamento all'istante i-1 da inserire nella (6 o nella (8. 

Se si specializzano la (2 e la (4 per t=0, ricavando velocitd ed accelerazione al tempo t=0, 
ricavando dalla (2 la x iU , e sostituendola nella (4, e ricavando da questa x,.,, si ottiene: 



x i-i = x i ~ AtXj - (Atf x t 



/2 



(9 
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E chiaro che le grandezze all'istante t=0 sono un dato del problema, in generale sono nulle. 



Vi sono due requisiti essenziali cui i metodi di integrazione al passo devono soddisfare: devono 
essere "stabili", e devono assicurare la "convergenza" alia soluzione esatta. II primo requisito e 
dovuto alia caratteristica che alcuni sistemi di integrazione numerica hanno di divergere 
indefinitamente se il passo di integrazione non e sufficientemente piccolo rispetto al periodo 
proprio della struttura. Piu precisamente se si calcola la risposta dell 'oscillatore semplice non 
smorzato nel caso di oscillazioni libere il cui risultato esatto e: Acos cot, nel caso del metodo delle 
differenze centrali la risposta pud risultare non oscillatoria ed anzi cresce indefinitamente se 
At>T/n. II metodo da comunque luogo a risposte sinusoidali ma con periodi inferiori rispetto a 
quello reale, questa tipo di errore, detto di convergenza, si riduce al diminuire del rapporto At/T 
come si vede in figura. 



o 

-0,05 
t -0,1 
b -0,15 
-0,2 

-0,25 J 



) "~~^5 


t\ 0 


2 0 



























dt/T 



Metodo delle differenze centrali, scostamento tra periodo proprio calcolato ed effettivo per 
oscillazioni libere non smorzate 



In pratica nel caso sismico la necessitd di rappresentare Vazione in modo accurato impone in 
genere di scegliere un passo di integrazione non superior e a 0.02 secondi, mentre valori tipici 
vanno da 0. 005 a 0. 01 secondi. Pertanto, nel caso di strutture civili, ad un grado di libertd, nelle 
quali il periodo proprio e in genere superior e a 0.1 secondi, il rapporto At/T e inferior e a 10, il 
metodo risulta pertanto stabile (At/T<l/n: condizione di stabilita) ed accurato (si veda la figura). 
Se la struttura entra in campo plastico i due requisiti divengono piit facili da soddisfare, in quanto 
I 'effetto globale e equivalente ad un allungamento del periodo. 



In definitiva i passi da seguire nella programmazione del metodo delle differenze centrali sono i 
seguenti: 

passo iniziale: valutazione dell'accelerazione al tempo t=0 (eq.(l, eq.(9, (8: 



x(t = 0) 



(-dx(t = 0) - kx(t = 0)-ma(t = 0)) 



m 
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m d 
[{Atf 2At_ 


Xj = -ma(t 0 ) - 


m d 
(Atf 2At_ 


2m 

X - 1 [ (At) 2 \ 






A tutti i passi successivi (eq.(8): 


m d 
[(At) 2 ' 2At_ 


x i+1 =-ma{t t )- 


m d 
[Atf 2 At 


j 2m 

r I w 


X i 


(11 


Se si devono calcolare le altre grandezze dipasso si utilizzano la (2 e la (4: 


■ _ X i+\ X i-l 
X- — 

2At 

x i+l - 2x i + x M 


(12 
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1.5.3 Metodi impliciti: il metodo di Newmark 

Si ipotizzi che si possa approssimare V accelerazione di risposta (relativa) in ogni passo con una 
accelerazione costante pari alia media delle accelerazioni di inizio (punto i) e fine passo (punto 
i+l). 

In tale ipotesi si possono ricavare le grandezze del moto lungo il passo ed a fine passo: 



.. x M +x t 
' " 2 


(13 




(14 


x m =x t +x t At+ \x M +x t ) 


(15 


La (14 pud essere scritta come 


x M = i,. + [(1 - 0.5)Af]x,. + (0.5At)x M 


(16 


x M = x i + x,At + [(0.5 - 0.25)A^ 2 Jx, + [o.25A? 2 Jx !+1 


(17 


Le espressioni (14a ed (14b possono essere generalizzate secondo le due espressioni seguenti 
dovute a Newmark (1): 


x m = x , + [(! - a)At^c, + {aAt)x M 


(18 


x M = x i + x t At + [(0.5 - P)At 2 \c i + \pAt 2 Jx ( . +1 


(19 



Se i valori di a e ft sono rispettivamente 0.5 e 0.25 le due espressioni coincidono con le (16 e (17. 
Si trova immediatamente che se i valori di a e (5 sono rispettivamente 0.5 e 1/6, le due espressioni 
coincidono con la soluzione relativa al caso in cui Vintegrazione nel passo si esegua utilizzando 
una variazione lineare della accelerazione relativa tra inizio e fine passo (secondo la ben nota 
regola dei trapezi). 

La (19 pud essere utilizzata per ricavare lo spostamento di fine passo. La (18 per ricavare la 
velocitd di fine passo. Si deve utilizzare un criterio iterativo assegnando inizialmente un valor e di 
tentativo alia accelerazione. Se a fine passo, sostituendo i valori trovati nella equazione di 
equilibrio dinamico (1 si ottiene un valor e della accelerazione sufficientemente prossimo a quello 
di tentativo, si pud ritenere di aver raggiunto la convergenza, altrimenti si deve iterare sino a 
convergenza. 

1.5.4 Riferimenti bibliografici 

Newmark,N.M., 1959, A Method of Computation for Structural Dynamics ", Journal of Engineering 
Mechanics Division, ASCE 85 
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1.6 Oscillatore Semplice: Stima delta risposta massima non lineare 

E' stato dimostrato che lo spostamento massimo di un oscillatore semplice elasto plastico soggetto 
ad un accelerogramma, se il periodo proprio e elevato, e circa uguale a quello dell 'oscillatore 
indefinitamente elastico. Pertanto il taglio massimo nella struttura e pari al taglio elastico diviso 
per la duttilitd massima richiesta. 

Se il periodo proprio e breve, la risposta in spostamento e tale da aver la stessa energia di 
deformazione massima dell 'oscillatore elastico indefinito. Detta E e i=l/2 Kx e , detta E ep =l/2 

2 2 

Kx y +1K (x u -x y ) = Kx tl Xy-l/2 Kxy , uguagUando i due termini si ha: 

x u =0.5(x e 2 / x y + x y )= 0.5(q x v + x v ) 

esprimendo x u normalizzato rispetto adx y si ha: 

H=x u /x y = 0.5[(x e /x y ) 2 + 1)]= 0.5(q 2 + 1) 
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In genere se Vazione sismica e rappresentata dallo spettro di risposta elastico di progetto, ed il 
periodo proprio si trova nella zona a velocitd costante: T>T C , si assume valida Vipotesi di ugul 
spostamento, al di sotto di tale periodo Vipotesi di ugual energia e ragionevolmente cautelativa. 
Detto q il rapporto tra forza di inerzia massima F e se la struttura rimane in campo elastico e la 
resistenza della struttura F y , qualora il periodo della struttura T* sia inferiore a T c si pud utilizzare 
V espressione : 



dT = 

max 



d 



e,max 



q 



z 



c 



r 



>d 



e.max 



per il calcolo della riposta massima in spostamento, nella espressione data d e>max e lo spostamento 
nel caso di risposta elastica. 
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d(q,0.8-Tc) 
d(q,0.7-Tc) 
d(q,0.6Tc) 
d(q,0.4-Tc) 
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2 Strutture a Piu Gradi di Liberia: Coordinate Generaiizzate 

Si consideri il telaio a 2 piani in figura. 




XI 



Si trascura la deformabilitd assiale e per taglio, le travi sono infinitamente rigide 
(Unitd dimisura: [F]=[kN], [L]=[m], [m]=[ton]) 



Caratteristiche geometriche 
Dimensione pilastri: 
Altezza di piano: 
Lunghezza campata: 
Numero di pilastri: 
Interasse telai: 

Momento d'inerzia dei pilastri: 

Rigidezza dei pilastri: 

Rigidezza di piano: 
Caratteristiche meccaniche 
Modulo elastico del els: 
Masse 

Massa impalcato: 

Massa sulla trove di piano: 

Massa di piano: 



bxh 
h p 
L 
n 



I=—bh 3 
p 12 



PP 



K 



kp n xkpp 



mi 

m t =miXi n 
m=m t xL 



=300x300 mm z 
=3,000 m 
=6,000 m 
=2 

=5,000 m 
=6,75 xlO' 4 ' m 

=9xl0 3 kNm 1 
=18xl0 3 kNm 1 

=30000 MPa (=30000x1 0 3 kNm 2 ) 

=0,70 ton m 2 
=3,5 ton m 1 
=21 ton 
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2. 1 EQUAZIONI DEL MOTO E MA TRICI DELLE MASSE E DELLE RIGIDEZZE 




■* *• 



77 sistema di equazioni che governano il moto della struttura si possono ricavare utilizzando 
differenti metodi : 

1) Equilibrio alia d'Alambert (scrittura diretta delle equazioni di equilibrio dinamico) 

2) Metodo energetico (Principio di Hamilton o della conservazione dell'energia totale del 
sistema) 

3) Elementi finiti 

Mentre il metodo 3) e utilizzato ampiamente nei programmi di calcolo perche si presta ad un 
elevato grado di automatizzazione, i primi due metodi sono il frutto del tradizionale approccio alia 
dinamica Lagrangiana. 

Con riferimento al telaio dell'esempio viene ora descritto il metodo 1). 

Nella figura accanto sono illustrate le sollecitazioni di taglio ad ogni piano (verso ed ampiezza). II 
telaio e soggetto ad uno spostamento alia base x g . Effettuando V equilibrio alia traslazione delle 
masse di piano possiamo scrivere il sistema di equazioni seguenti, ognuna delle quali mostra come 
la forza d'inerzia (proporzionale all'accelerazione assoluta della massa) equilibri la reazione 
elastica (proporzionale alia rigidezza di piano kp e dipendente dallo spostamento relativo di 
piano). 



Jm(Xj +x g )- k p (x 2 - x l ) + k p x { = 0 

[m(x 2 +x g ) + k p (x 2 -x l ) = 0 ^ 
che si pud riscrivere mettendo a fattor comune le incognite xj e X2 ossia gli spostamento relativi 
piano (vedi figura). 
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mx l + 2k p x l - k p x 2 = -mx 



mx 2 - k x x + k p x 2 



-mx„ 



Volendo scrivere le precedenti in forma compatta si pud utilizzare la notazione matriciale: 



m 0 
0 m 



Ik.. 



m 0 
0 m 




MX+DX+KX = -MIx g 

nella quale M e la matrice delle masse, K e la matrice delle rigidezze, I e il vettore di 
trascinamento edX e il vettore incognito degli spostamenti di piano. II termine a secondo membro 
rappresenta il vettore delle forze equivalenti al sisma. 

2.2 SOLUZIONE CON APPLICAZIONE AD UN TELAIO PIANO 

Nel caso dell'esempio le matrici delle masse e delle rigidezze assumono le espressioni: 

~m 0] [21 0" 

XI - 

0 m 0 21 



XI 



matrice delle masse 



K = 



2k, 



K = 



3.6 -10 4 
-1.8-10 4 



-1.8 10 4 
1.8-10 4 



matrice delle rigidezze 



Se si ipotizza di conoscere una soluzione del tipo X=</>y(t), ovey(t) e uno scalar e che viene detto la 
coordinata generalizzata e serve a modulare il vettore </>, quest 'ultimo da la forma della 
deformazione. 
Ad esempio: 
"0.115" 



0.186 




Spostamenti del primo e secondo piano 
Per il principio dei lavori virtuali il lavoro delle forze interne ed esterne fatto in uno spostamento 
virtuale deve essere nullo perche vi sia equilibrio. 

Sostituendo a X V espressione dello spostamento X=0y e le sue derivate temporali, assunto proprio 
oX=<f> dy quale spostamento virtuale, premoltiplicando per $ i termini dell' equazione di equilibrio 
si ottiene: 

dy ■ ftMfy + dy ■ fM\x g + dy • fD^y + dy ■ 0 T K0y = 0 

Dividendo per dy si ottiene: 

<f) T M<f)y + <f> T MIx + <f) r D</>y + <fi T K<f>y = 0 
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Che pud riscriversi come: 

<f> T M<j)y + 0 T D<f>y + 0 r K<f>y = -<f> T Mlx g 

Da cui: 

(j/Mfy + <f> T D<f>y + fK<f>y = -0 T Mtt g 

Ponendo: 
4> T M(I) = m 

fD0 = d 

fK0 = k 

fMl = L 

L 'equazione di equilibrio si riscrive: 
my + dy + ky = -Lx g 

m = massa generalizzata, d = smorzamento generalizzato, k = rigidezza generalizzata. 
Dividendo per la massa generalizzata si ottiene: 
y + Ivay + a> 2 y = -px g 

Ove p=L/m, e detto coefficiente di partecipazione. 

L 'equazione e quella di un oscillatore semplice soggetto ad una accelerazione amplificata p volte 
rispetto a quella effettiva. 
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3 Strutture a Piu Gradi di Liberia: Analisi Modale 



3. 1 CALCOLO DELLE FREQUENZE E PERIODI PROPRI DEL SISTEMA 



Si consideri ora il caso nel quale sia assente la forzante esterna in genere denominato problema 
delle OSCILLAZIONI LIBERE 
MX + KX = 0 

Si pud ipotizzare una soluzione di tentativo del tipo 
X := (|)-sincot 

Sostituendo la precedente nel sistema di equazioni si ottiene un nuovo sistema nel quale le 
incognite sono il vettore delle ampiezze <j) e la frequenza del moto co, la determinazione della cui 
soluzione e delta PROBLEMA AGLI AUTOVALORI 

(K-co 2 M)(j) = 0 

E4> = 0 

ossia in forma compatta 

II sistema cosi ottenuto e un sistema omogeneo che come ben noto dalla matematica ammette 
soluzione (autosoluzione) se e soltanto se il determinante della matrice dei coefficienti possiede 
determinante nullo. Dunque per calcolare le autosoluzioni del sistema occorre imporre la 
condizione: 
det[E] = 0 

Nel caso del telaio in esame il problema agli autovalori si scrive semplicemente: 



^ 2kp kp ^ 
V _k P k P J 



- co 



m 0 
0 m 



<t>12 



ossia 

f 

2k n 



2 

co m 



" k P 



" k P 



CO 



<t>12 



Imponendo Vannullamento del determinante della matrice E si ottiene la seguente equazione 
algebrica biquadratica nell'incognita co 
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( 2 V 2 ^ 2 
I 2kp - co m J\ kp - co m ) - k p := 0 

2 4 2 2 

m co - 3mkpC0 + kp := 0 

le cui soluzioni sono date dalle espressioni seguenti 
2 3mkp + 



coi 



C02 



9m 2 k p 2 



2 2 
4m k p 



2m 



3mkp - 



9m 2 k p 2 



2 2 
4m k p 



2m 



/ due valori col e a>2 sono detti FREQUENZE PROPRIE del sistema 
Calcolo coefficient 

Polinomio caratteristico ax^+bx+c=0: 



a := m 



b := -(3-kn-m) 

aaa VP/ 



c :=k n 

AAA V 



-b - 



b -4-a-c 



2- a 



Xl = 327.399 

b + yi 



x 2 : 



4-a-c 



2- a 



X 2 = 2.244 x 10 
FREQUENZE PROPRIE 
PERIODI PROPRI 



coi 




coi 


= 18.094 




2-71 


Ti 






K>1 


Tl 


= 0.347 


co 2 


= fil 


C02 


= 47.371 




2-71 


T 2 






co 2 


T 2 


= 0.133 



- Uso del problema agli autovalori generalizzato implementato in Mathcad 
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Spesso i programmi di calcolo automatico sfruttano metodi numerici approssimati per la 
valutazione della frequenze proprie 
FQUAD := genvals (K, M) 



f 



FQUAD = 



2.244x 10' 
v 327.399 , 



O :=4 FQUAD 
T := 2-;t Q~ 1 

AAA/ 

n.f 4 " 71 ] 

^18.094^ 
(FREQUENZE PROPRIE) 
'0.133 N 

^0.347 y 
(PERIODI PROPRI) 
CALCOLO DEGLI AUTOVETTORI 
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Poiche gli autovettori sono definiti a meno di una costante e possibile scegliere differenti condizioni 

di normalizzazione 

1 ) metodo della componente unitaria 

Imponendo la condizione che almeno una componente di ciascun autovettore sia 1 e partizionando 

la matrice del problema agli autovalori si ha 

EH*:=0 

problema agli autovalori 
autovettore i-mo 

Ej := K - Xi-M 

matrice del problema agli autovalori 
partizionamento della matrice del problema agli autovalori 
f en Ei,l(A 

^,01 Ei,ll ) 

se si pone ad esempio 
fcl := 1 

dalla seconda equazione si ricava il valore di <j)l2 

poiche gli autovalori hanno reso dipendenti le due equazioni la prima sard automaticamente 
soddisfatta. Nel caso di piu di due gradi di libertd il termine E;,01 rappresenta la prima colonna 

della matrice dei coefficienti moltiplicato per 1, primo valore dell' autovettore posto =1, esso fa si 
che i termini noti del problema siano proprio i valori e,,01 

passiamo al nostro esempio 

Autovettore 0J 
4>n:=l 

Ei :=K-Xi-M 
Eli:=El M 

<t>21 — E l 1 •% 



<l>la : = 
<t>la = 



l 

1.618 y 
Autovettore <j>2 
<I>12 : = 1 
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E 2 := K - X 2 -M 
E22 := E2 



1 , 1 



E02 := E2 Q x 
<t>22 := ~ E 22 



E02 



tea 



tea : 



J 



<t>22 
-0.618 



La forma generate prima ricavata pud essere particolarizzata al caso attuale. Ad esempio con 
riferimento al primo modo la prima equazione del moto si pud scrivere: 
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(2-kp - X r m)-bu ~ V<l>12 == 0 
dalla quale si evince, che posto 

2-kp - X\-m 
(t> 12 = 1.618 

Analogamente per Vautovettore del secondo modo 

2 ) metodo della normalizzazione della matrice delle masse 

In tale caso si impongono le seguenti condizioni: 

*i T M^:=l 
i:= l..n 
n=n° piani 

Raccogliendo gli autovettori <jn in un 'unica matrice 0 detta matrice modale le precedenti possono 
scriversi in maniera compatta 

d> T M O :=] 

dove la matrice I e la matrice unitd 
1 0 s 

I: 



.0 I, 

Per il caso dell'esempio imponendo la precedente al problema agli autovalori prima definito si 
ottengono le seguenti componenti del i-mo autovettore 



2 



(2kp - Xi-mj + kp 



2 



te- 
rn 

Calcoliamoci dunque le componenti dei singoli autovettori 
Autovettore <f)l 



+1 = 



Xp 2 



(2kp - Xymj + kp^ 



2 

1 -m$n 



m 



V *12y 
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Autovettore <f)2 




m 



/ \Z 

(2k p - l 2 -mj + k p 



2 



1 - m-(|>21 



m 



<l>2 : 



1*22 J 
f 0.186 



to 



-0.115 



NB. Una maniera alternativa di imporre la normalizzazione rispetto alia matrice delle masse e 
quella di utilizzare i vettori normalizzati col metodo della componente unitaria e poi sfruttare la 
proprieta che gli autovettori sono definiti a meno di una costante per valutare i nuovi vettori. In 
particolare ponendo le seguenti condizioni 
to := Crtoa 

to := Q'toa 
T 

(|>1 -M (|>i := 1 
T 

§2 -M-(|>2 := 1 

(normalizzazione rispetto alia matrice delle masse) 

dove i vettori 

toa 

toa 

sono i vettori ricavati col metodo della componente unitaria 
si ottengono le soluzioni per i coefficienti CI e C2 



1 


toa T 


M-toa 


1 


toa T 


M-toa 



/ nuovi vettori <f)l e <f>2 risultano ora normalizzati rispetto alia matrice della masse 
Nel caso dell'esempio si ha: 
Ci = 0.115 
C 2 = 0.186 
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<t>l := Cr<|>i a 





0 


1 


0 


1 




1 


2 







( 0.1 15^ 






v 0.186j 




( 0.186 


4>2 = 




V -0.115 J 



che sono esattamente i vettori precedentemente determinati imponendo direttamente la 
normalizzazione rispetto alia matrice delle masse. 

Naturalmente per valutare gli autovettori e possibile, come logico, utilizzare anche un comando 
automatico come quello implementato in Matchad 

Mathcad usa una normalizzazione differente dalle due esposte precedentemente 
O := genvecs (M,K) 



'0.526 0.851 ^ 
v 0.851 -0.526J 
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-0.2 o 0.2 



1° modo di vibrare 
11° modo di vibrare 
Ti = 0.347 

T 2 = 0.133 

NB : Come e usuale il primo modo di vibrare risulta lineare mentre il secondo risulta intrecciato 
con valore negativo dello spostamento in sommitd 
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3.1.1 PROPRIETA ' DI ORTONORMALITA ' DEGLI AUTOVETTORI 

A perscindere dal metodo di normalizzazione adottato si dimostra facilmente che gli autovettori 
appena ricavati possiedono una importantissima proprietd detta di ortonormalita che si esprime 
come segue 

se 

i e diverso da j 

fc T .K-4>j:=C 
se 

i e diverso da j 

Se poi si adottasse la normalizzazione degli autovettori rispetto alia matrice delle masse si 
otterrebbero le condizioni aggiuntive seguenti 

*i T -M-4i:=l 

T 2 
fa -K(|)j:=coi 

Ossia in forma compatta: 

<J> T -M-<1> :=] 

f 1 0 (A 



I := 



0 .. 0 
0 \) 



O K O := A 

f 2 
t»i 



A := 



0 



0 co 



n / 



Con le precedenti e possibile rendere le equazioni del sistema indipendenti tra loro (vedi calcolo 
della risposta col metodo della sovrapposizione modale). 

Nel caso dell'esempio la verifica dell' ortonormalita degli autovettori e immediata : 



§\ M-<|>2 = 0 
0 



T 

<)>! M-(|>2 



T 

<P2 K '<t>l : 



-2.842 x 10 



14 



13 



2.274x 10 

Calcolo a posteriori delle frequenze proprie del sistema 
A:= fo T M-4>i 



B:=(|)2 M-tyi 
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C:=<bi K-cbi 



D:=*2 K-<|i2 

- 

(Oi = 18.094 
Frequenza 1° Modo 

co? := — 

(WM | g 

©2 = 47.371 

Frequenza 11° Modo 
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3. 1.2 CALCOLO ITERATIVO DEI MODI E FREQUENZE PROPRIE 

Si illustra un metodo iterativo per il calcolo di autovalori ed autovettori, utilizzabile con qualsiasi 
programma statico, utile per il controllo dei risultati dell'analisi dinamica. 
data I'espressione agli autovalori: 
(K-co 2 M)4) = 0 
se si premoltiplica per 

-\k- 1 (K-6) 2 M)0 = O 

CO 

si ottiene I'espressione: 
(\l-K- 1 MU = 0 

CO 

che da luogo alia forma iterativa: 

CO 

0n l =K- 1 M^°;^n 1 =Xr$ 

CO 

o anche 

ove con I'apice 0 ed 1 si indica iterazione 0 e iterazione 1. Ciod si assegna una forma di tentativo e 
si ricava dall'espressione detta una forma successiva. Poiche la frequenza e incognita essa pud 
essere inclusa nel vettore § che deve essere poi normalizzato. La normalizzazione si pud fare a 

posteriori osservando che il rapporto tra il valore effettivo di un autovettore in un punto e quello ed 

il suo valore normalizzato e appunto co . Pertanto eseguite un certo numero di iterazioni e possibile 
calcolare il quadrato della frequenza come il valore prima dell'iterazione e dopo I'iterazione ma da 
normalizzare: 

co 

da questa espressione si ha che il valore della frequenza si pud ottenere come rapporto tra i valori 
che I'autovettore alia iterazione precedente e quello successivo hanno in un punto preciso della 
struttura. Da punto a punto tale rapporto varia, il valore esatto sard compreso tra il massimo ed il 
minimo dei rapporti trovati, ovvero, pesando i valori con la distribuzione delle masse: 

2 d>n l 'Md>° 
co =— — 

(j>n M<j)n 

il metodo converge rapidamente, ci siferma quando a iterazioni successive i valori ottenuti sono 
simili. 

ad esempio nel caso in questione si ha, con I'assunzione di un autovettore iniziale unitario ai due 
piani: 
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(<|ni):=K l -M-h 

( 

2.333x 10 ' 

fai = 

y 3.5 x 10 j 
normalizzando a 1 la prima componente si ha: 
¥\ 

¥\ ■= 



i 



f -3^ 

1.944x 10 



col 



^3.111x 10 J 



,0.5 



V 



T 



col = 18.096 

Tl :=2-7t— 
col 

Tl = 0.347 

v«o/ rendere il vettore normalizzato in modo ortonormale: 

T 

c 

0.116 
0.185 



(|>1 M(|)i = 1 

3.2 CALCOLO DELLA RISPOSTA: METODO DELLA SOVRAPPOSIZIONE MODALE 

Equazioni modali e fattori di partecipazione 

Una volta valutati i modi di vibrare e le frequenze proprie del sistema e possibile calcolare la 
risposta sismica del nostro telaio. Spesso in fase progettuale si adotta quello che viene usualmente 
denominato metodo della sovrapposizione modale. 

L'idea base del metodo e legata alia proprieta di ortonormalita degli autovettori e alia 
normalizzazione degli stessi rispetto alia matrice delle masse che rendono disaccopiate le equazioni 
del moto. Infatti siponga: 
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X:=(D\ 

dove Y e il vettore delle coordinate normali mentre X e il vettore degli spostamenti di piano e 0 e 
la matrice modale 
$ :=(<h te <M 

essendo n il numero dei gradi di libertd 

Sostituendo le precedenti nelle equazioni del moto e premoltiplicando ambo i membri per il vettore 
<f)i trasposto si ottiene un nuovo sistema in Y, composto da n equazioni disaccopiate, la cui generica 
espressione e la seguente: 

O t MOY + O t KOY = -O x MIx g 

i:= l..n 

Con la normalizzazione dei modi rispetto alia matrice delle masse l'i-ma equazione diventa 
Yi + C0i 2 Yi = -Pi^g 



ossia V equazione di un generico oscillator e semplice di frequenza coi soggetto ad una accelerazione 
moltiplicata per il fattore pf. Quest'ultimo e detto fattore di partecipazione ed e definito come segue 

pi := eh -M-I 

e rappresenta il lavoro compiuto dalle forze sismiche 
MI 

nel modo 
<h 

E' allora chiaro che se pi e alto c'e un forte coinvolgimento del modo <jA nel moto imposto dal 
sisma. E' altrettanto chiaro che piu il modo e intrecciato e piii c'e la possibilitd che il relativo 
fattore di partecipazione sia piccolo. 
Smorzamento Strutturale 

E' usuale considerare il singolo oscillatore smorzato con fattore di smorzamento dp. In tal caso 
I'equazione diventa 

In genere il valore dello smorzamento modale viene determinato immaginando una matrice degli 
smorzamenti C cosi' definita 
C:= a-M + b-K 

AAA/ 

Smorzamento alia Rayleigh 

I parametri a e b vengono definiti sulla base dei due modi piu importanti: 

T T T 

(|>i C-([)i := a <)>i -M-t^i + b-<(>i -K-^j 

T T T 

(|>j -C-q>j := a- 4>j -M-ty + b - cf>j -K-(|>j 

Ponendo 

T 

(ft -C-(|>i := 2-^i-cOj 
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si ottengono le espressioni dei coefficienti a e b 
a:=2.& 



(X»i + COj 



b:=iy 

COj + 03[ 

Massa eccitata 

NB : la somma dei quadrati dei fattori di partecipazione e, per effetto delle ortonormalizzazione 
degli autovettori, uguale alia massa totale del sistema: 

Pl 2+ P2 2=m tot 

Valutata la risposta del singolo oscillatore e suffciente tornare al vettore X per determinare la 
risposta in coordinate lagrangiane. Lo spostamento di piano pud essere espresso semplicemente 
come somma dei contributi dei singoli modi (sovrapposizione modale) 

n 

j=i 

dove n e il numero dei modi 

In genere e possibile considerare nella sommatoria solo i modi di vibrare piii significativi. A tale 
scopo e utile definire la massa eccitata dal singolo modo: 

2 

Pi 

e; := 



mtot 

essendo mtot la massa totale dell'edificio 

Valori bassi sono indice di un modo poco eccitato e che contribiusce poco alia risposta ( valore 
basso nel termine della sommatoria) . Dunque la massa eccitata e un mezzo per avere un'immediata 
idea di quanti modi considerare per avere una risposta attendibile. 
Nel caso dell'esempio si ha: 

vettore delle incidenze per moto orizzontale, detto anche vettore di trascinamento: 

<) 

Fattori di partecipazione 

Pi := (fr^M-l 
Pi =6.315 

Fattore di partecipazione 1° Modo di Vibrare 

T 

P2 := 4>2 M-I 
p 2 = 1.489 

Fattore di partecipazione IF Modo di Vibrare 

PT:=P1 + P2 

pT = 42.093 
77 



Prof. Camillo Nuti Dispense di Costruzioni in Zona Sismica 2007 



2-m = 42 

si noti il basso valoredel secondo fattore di partecipazione sintomo del fatto che il secondo modo e 
relativamente (al primo modo) poco attivato. A tale scopo e utilie anche il calcolo della massa 
eccitata dei singoli modi di vibrare. 

2 

PI 



£1 


:= IOC 

2-m 


£ 1 


= 94.944 


% 






2 

P2 

:= IOC 

2-m 


^2 


£2 


= 5.279 


% 



La risposta del singolo oscillatore pud essere determinata in vari modi ad esempio 
1 ) Calcolo della risposta istante per istante con llntegrale di Duhamel 
2 ) Calcolo della risposta massima col Metodo dello spettro di risposta 
Metodo dell' Integr ale di Duhamel 

x . = 2> Jn-Pi*! — e ^ ,( " T) sino^t-^dt 
Metodo dello Spettro di Risposta 

Se in presenza di una accelerazione ag lo spettro di risposta in spostamento e pari 
Sd(^i^i) 

la risposta massima del singolo oscillatore modale sottopostoi ad un accelerazione pi* ag sard 

yima X :=Pi-Sd(^i^i) 

Passando ora in coordinate non normali il vettore degli spostamenti di piano pud essere 
determinato con un opportuna combinazione che in genere viene scelta come segue 

Xmax := ^^rtmy 

dove il vettore Ymax racchiude i valori massimi delle risposte dei singoli oscillatori modali 
Sicche l'i-ma componente di Xvarra 

Ximax^ / ^ ifyima? 

Usualmente la scelta del numero di modi da considerare dipende dalla % di massa eccitata dei 
singoli modi. 
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3.2. 1 ESEMPIO: Calcolo della risposta di un telaio con lo spettro di risposta 

Spettro di progetto nuova normativa: zona 1 



a g max:= 0.339.81 
a gmax = 3.434 



Caratteristiche dello spettro di risposta elastico (Categoria suolo A) 



S:=l 

AAA 

T B :=0.1f 
T C := 0.4C 
T D :=2.C 
r, :=1 

(Smorzamento 5%) 

q:=4 

(Fattore di Struttura) 

Spettro di risposta elastico in accelerazione 



S a (T) := 



gmax 



1 + (ri-2.5- l) 



if 0 < T < Tg 



gmax 



Sri-2.5 if T B <T<T C 



gmax 



S-r|-2.5- 



if T C < T < T D 



V i J 



gmax 



S-T|-2.5' 



T C Tp 

V T 2 



if T D < T 



S a (0.4) =2.146 
S a (l) =0.858 

Spettro di risposta in spostamento 



S d (T):=^— J -S a (T) 

Nel caso del telaio dell'esempio gli spettri in spostamento e accelerazione dei singoli oscillatori 
modali varranno: 

Sdl :=Pl-Sd(Ti) 

S d i = 0.041 

Spostamento massimo oscillatore modale 1 

S a (Tl) 

S a l : = Pi 

9.81 

S a i = 1-381 

Accelerazione massima oscillatore modale 1 (in g) 
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Sd2 :=P2'Sd(T2) 
S d2 = 1.325 x 10~ 3 

Spostamento massimo oscillatore modale 2 

S a(T 2 ) 

Sa2:=P2 ^r 

S a2 = 0.303 

Accelerazione massima oscillatore modale 2 (in g) 
E dunque la risposta massima ai singoli piani varrd 
Vettore degli spostamenti di piano 
( 0.041 

Y:= 
O -- 



V 0.001325 
r 0.116 0.186 



X:=0\ 



0.185 -0.115 



X = 



f -3^ 

5.033x 10 
v 7.507x 10" 3 j 



In modo formalmente equivalente si pud scrivere: 

■®>A x ^ « 

Xmmax:= > u 
n =0 

Ove con u(i) si sono indicati i vettori colonna dei contributi dei singoli modi. 



3.2.2 Criteri di sovrapposizione modale con I'uso dello spettro di risposta 

La formulazione sopra riportata e corretta se applicata alia risposta Y(t), e una stima per eccesso 
della risposta, infatti i massimi delle due coordinate Ynon sono contemporanei. 
Si adotta usualmente uno dei due metodi di seguito indicati: 

E = (S Er) 12 

E = (I^j pjj Ei Ej) 1 * 

Dove con E e indicato il contributo di ciascun modo alia risposta. 

il primo criterio si applica come riportato di seguito: 
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Xmax:= 



n =0 



0.5 



Xmax= 



4.793x 10 



7.661 x 10 



// secondo criterio di sovrapposizione modale e piu preciso, tiene conto della correlazione tra i 
modi, dovuta alia vicinanza delle frequenze ed alio smorzamento. La prima formulazione della 
correlazione modale e dovuta a Rosemlueth ed Elourdy, cosi come la formula di combinazione. Qui 
si mostra la formulazione suggerita da Der Kiureghian : 



P :=■ 



H (i + P)-P 



p = 8.856x 10 



(l - P 2 )\ 4.^.(1 + pf 



Xmaxl: 



i (u«r 



n =0 



+ p -uv 



0.5 



C 



Xmax 1 - 



-3 



4.794x 10 



V 7.66x 10 J 



Mettiamo a confronto i tre risultati ottenuti con i tre criteri di sovrapposizione modale: 



f 



X 



5.033x 10 



7.507x 10 



Xmax= 



4.793 x 10 



7.661 x 10 



Xmaxl = 



4.794x 10 



7.66x 10 



Si osseva come il primo criterio, eccessivamente cautelativo, dia risultati non particolarmente piu 
gravosi in questo caso. Esso rappresenta una stima per eccesso. II secondo ed il terzo criterio 
danno luogo a risultati pressoche coincidenti. Questo fatto e dovuto alia scarsa correlazione 
modale. Infatti i periodi propri sono piuttosto lontani: 0.35 e 0.14 secondi. La correlazione e 
sempre pressoche nulla al di sopra di una distanza di 0. 1 secondi. 
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Correlazione modale al variare del rapporto tra i periodi: f3=Ti/T2 
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3.3 TELAIO SHEAR-TYPE 2 PIANI ISOLA TO ALLA BASE 





PREDIMENSIONAMENTO SISTEMA DI ISOLAMENTO 
periodo diprogetto struttura isolata 



Tiso :— 2 



Massa base isolamento 
m iso := m 

miso = 21 

Massa totale del sistema di isolamento e della sovrastruttura 
mtot :=2 m+ m iso 
mtot = 63 

Calcolo rigidezza della base isolata considerando la sovrastruttura rigida 



4-n 



-mtot 



k iso = 621.785 



3.3. 1 MATRICI DEL SISTEMA DI EQUAZIONI DEL MOTO 



M := 

/WW 



M 



K:= 



m 



0 0^ 



0 m 0 
0 0 mj 

(21 o o^ 

0 21 0 

0 0 21 
(kp + kj so 



2k n 



0 ^ 
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K 



1.862x 10 -1.8 x 10 0 
-1.8 x 10 4 3.6 x 10 4 -1.8 x 10 4 



V 



-1.8x 10 1.8x 10 J 



J 



3.3.2 CALCOLO DELLE FREQUENZE E PERIOD! PROPRI DEL TELAIO 

- Uso del problema agli autovalori generalizzato implementato in Mathcad 
FQUAD := genvals (K, M) 

AAAAAAAAAAAA/ ° \ 7 ■/ 

( 3^ 

2.576x 10 
FQUAD = 872073 

v 9.682 
(FREQUENZE PROPRIE) 
Q :=Vfquad 
( 50.759^1 

29.531 

^3.112; 

T := 2-7t-Q~ 1 

AAA/ 

C 0.124^ 
0.213 
2.019; 



(PERIODI PROPRI) 

Si noti come la prima frequenza sia prossima alia frequenza del sistema considerando la 
sovrastruttura rigida e che la seconda frequenza sia relativamente piu grande 



ISO 

K iso 
3m 



coi so = 3.142 



Wis 



o>bas '■= 

m 

«bas = 5.441 
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Nel caso di azione sismica questo fatto risulta estremamente rilevante in quanto I'aumento del 
periodo proprio della struttura porta la struttura verso zone dello spettro di progetto con ordinate 
piii basse, con conseguenti minori azioni sulla struttura. 
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3.3.3 CALCOLO DEGLI AUTOVALORI 

1 ) metodo della componente unitaria 

Imponendo la condizione che almeno una componente di ciascun autovettore sia 1 si ha: 
Autovettore 0J 

e^=k-(q 2)0 ) 2 m 

]fofc= submatrix(Ei,l,2,l,2) 
Ejo := submatrix(Ei , 1 , 2, 0, 0) 

to=-En -E 10 
<t>ll 



hi 



<t>21 



0,0 

<t>21 

1,0 J 

1 > 



1.023 
Autovettore <f>2 

E iv :=K-(Q 1)0 ) 2 M 
E^- submatrix(Ei,l,2,l,2) 
E^:= submatri^E!, 1,2,0,0) 

<t>ll 




0.017 
V -0.983^ 
Autovettore <f>3 

^=K-(n 0)0 ) 2 -M 

E^- submatrix(Ei,l,2,l,2) 
E^:= submatri^E!, 1,2,0,0) 
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2 ) metodo della normalizzazione della matrice delle masse 



te 



' 0.1 24^ 
0.126 
v 0.128 y 



1 



tea -M-(|)2a 
^ :=C 2-*2a 



0.156 

2.665 x 10" 3 
-0.153 J 



1 



<l>2 = 
C 3 := , 

tea 1 -M- *l>3a 

te : = C 3-tea 
f 0.09 A 

(|)3= -0.178 

v 0.089 j 





0 


1 


0 


1 




1 


2 
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3 



i 2 



1 ' ' ' 

-0.2 0 0.2 




77° moJo (// vibrare 

NB : II primo modo e tale che la sovrastruttura praticamente si muova rigidamente sul basamento 
isolato e che tutta la deformazione sia concentrata negli isolatori. Come si dimostrerd nel 
successivo calcolo delle masse eccitate ad esso e associata la quasi totalitd della massa eccitata e 
dunque nella sovrastruttura giungono sollecitazioni assai ridotte rispetto al caso non isolato 



3 



i 2 




1 1 1 

-0.2 0 0.2 



fa 
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111° modo di vibrare 



3.3.4 CALCOLO FATTORI DI PARTECIPAZIONE 

vettore delle incidenze per moto orizzontale 



I := 

A/V 



1 



vettore deifattori di partecipazione 
pi =7.936 

P2 = 0.111 

P3 :=(|)3 T M-] 
p 3 = 0.022 

si noti il basso valore del secondo indice di partecipazione sintomo del fatto che il secondo modo e 
relativamente (al primo modo) poco attivato. A tale scopo e utilie anche il calcolo della massa 
eccitata dal modo i-mo. 

verifichiamo che la somma dei quadrati deifattori di partecipazione sia la massa totale del sistema 
isolato 

ryi 2 j 2^2 

J{&=Pl +P2 +P3 
pT = 63 

3.3.5 CALCOLO DELLA MASSE ECCITATE 

per definizione la massa eccitata del modo i-mo e il rapporto tra il quadrato del fattore di 
partecipazione e la massa totale del sistema 





^lioc 

3-m 


/wU/ 


ei = 99.98 


NB : 


la massa 


/v\2y 


2 

P2 

10C 

3-m 


£2 = 


0.02 


S3 : = 


2 

P3 

10C 

3-m 


S3 = 


7.523 x 10" 
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si noti come il primo modo rappresenti difatto il 100% della mass a eccitata, questo e abbastanza 
tipico dei sistemi ben isolati 

3.3.6 CALCOLO DELLA RISPOSTA SISMICA : Metodo dello Spettro di Risposta 



Spettro diprogetto nuova normativa: zona 1 
0.339 
3.434 



iW :=0 - 359 - 81 



a gmax~ 

Caratteristiche dello spettro di risposta elastico (Categoria suolo A) 
S := 1 

AAA 

TCJ= 0.4C 

(Smorzamento 5%) 

AAA* ^ 

(Fattore di Struttura: deve restare elastica) 
Spettro di risposta elastico in accelerazione 



gmax 



1 + Gl-2.5- l) 



if 0 < T < T B 



gmax 



Sri-2.5 if T B <T<T C 



gmax 



gmax 



S-ti-2.5 



S-ti-2.5- 



— 1 
I T J 



if T C < T < T D 



T C T D 



if T D < T 



V T j 

Spettro di risposta in spostamento 

v 2 



MP := [— J - S a(T) 

Nel caso del telaio dell'esempio gli spettri in spostamento e accelerazione dei singoli oscillatori 
modali varranno: 
,§dJU : =PrSd( T 0 ,o) 
S d i = 0.024 

Spostamento massimo oscillatore modale 1 

s , = Mv) 

w&kK Pi 9 81 
S a i = 6.216 

Accelerazione massima oscillatore modale 1 (in g) 
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to=P2-Sd(T lj0 ) 



S d2 = L092x 10 

Spostamento massimo oscillatore modale 2 



a (\o) 

51 



S a2 = 0.097 

Accelerazione massima oscillatore modale 2 (in g) 

Sd3 := P3- s d( T 2 5 o) 
S d3 = 3.787 x 10" 3 

Spostamento massimo oscillatore modale 2 

S a( T 2, 0 ) 

Sa3:=P3 ^r 

S a3 = 3.737 x 10~ 3 

Accelerazione massima oscillatore modale 2 (in g) 



E dunque la risposta massima ai singoli piani vale 





Sdl 


Y, := 


Sd2 


Y - 
/w2 " 


Sd2 






O := 




IWIA 





V 



0 0 0 



1 1 
(bo d)3 



0.024 > 
1.092x 10" 3 



1.092x 10 ) 
( 0.124 0.156 



0.09 



0.126 2.665 x 10 -0.178 
v 0.128 -0.153 0.089 ; 



Vettore degli spostamenti di piano 



X:=<D\ 

AAAA 
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X : 



3.193x 10 



2.802x 10 



^ 2.957 x 10 J 

In effetti in generate il sistema di isolamento da luogo ad uno smorzamento: <%i, superiore al 5%. Se 
si assume : ^i= 12%, che agisce sul solo modo, la risposta pud valutarsi riducendo la risposta 
spettrale moltiplicando i valori dello spettro di risposta, relativi ai modi interessati dalla 
distorsione del sistema di isolamento, per la seguente espressione: 



10 



12 

r[ = 0.767 

S^:=iVPrS d (T 0)0 ) 
Sdl =0.018 



9.81 



S a i = 4.768 

E dunque la risposta massima ai singoli piani, nel caso di smorzamento del sistema di isolamento 
pari al 12% vale 



Y„ 


= S dl 


Y, 


= S d2 


Y 
/w2 


= s d2 



0.018 

1.092x 10~ 3 

^1.092x 10" 3 J 
( 0.124 0.156 

0.126 2.665 x 10~ 
V0.128 -0.153 



0.09 



-0.178 
0.089 ; 



Vettore degli spostamenti di piano con smorzamento dell'isolamento=12% 
X:=<D\ 

AAAA 

f 



X : 



2.512x 10 



2.104x 10 



2.251x 10 J 
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Anche in questo caso valgono le considerazioni fatte per la sovrapposizione modale nel caso del 
telaio convenzionale. 

L'applicazione del metodo CQC (composizione quadratica completa) richiede I'uso della 
espressione della correlazione relativa al caso di smorzamenti differenti ( si veda il file Mathcad 
per i dettagli e le differenze di risultato. 

L 'espressione completa dovuta a Der Kiureghian e la seguente: 



8-(^o) ' 5 (§0+ |-P0l)-P01 



p01;= 

/VWVvV 



(l - P01 2 ) + 4-(S-Ef))-p01-(l + P01 2 ) + 4-($ 2 + £, 2 )-p01 2 



pOl = 0.069 
p02 = 5.426x 10~ 



Nella formula (301 e il rapporto tra la frequenza del primo e del secondo modo, £,0 e lo 
smorzamento del primo modo, dovuto essenzialmente agli isolatori, pari al 12%. 



r 



Xmmax= 



2.413x 10 



2.298 x 10 



V2.155x 10 J 



Xmax= 



f -3 > 

2.25 x 10 

2.296x 10~" 

^ 2.328 x W*) 



Xmaxl = 



f -3> 
2.28 lx 10 

2.392x 10~ 3 



V 



2.346x 10 



Si veda come in questo caso la correlazione modale sia leggermente piii influente per effetto dello 
smorzamento del 12% del primo modo. L 'andamento della correlazione modale p al cariare del 
rapporto tra i periodi e riportato nella figura in basso. Nella successiva e riportata la variazione di 
p se un modo ha smorzamento 5% al variare dello smorzamento dell'altro se il rapporto tra le 
frequenze modali vale 0.6 o il reciproco 1.67. Si vede che I'incremento dello smorzamento ha una 
influenza sensibile. Nel caso dell 'isolamento e possibile che i modi isolati abbiano frequenze simili 
se non addirittura coincidenti. In tal caso il valore di p pud essere elevato e la sovrapposizione di 
tipo CQC diviene particolarmente valida se non indispensabile. 
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0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

P 



Funzione di correlazione p al variare del rapporto f3 tra le frequenze modali: p- per 
smorzamento costante nei due modi pari a 0.05, pO- per smorzamento pari I 5% in un modo 
ed al 12% nell'altro. 



0.2 



0.15 



p01(^0) 0.1 



0.05 




0' 1 1 

0 0.05 0.1 0.15 0.2 

^0 

Funzione di correlazione p per un rapporto (3 tra le frequenze modali pari a 0.6 :, pO - per 
smorzamento pari al 5% in un modo ed al 12% nell'altro 



94 



Prof. Camillo Nuti Dispense di Costruzioni in Zona Sismica 2007 




95 



Prof. Camillo Nuti Dispense di Costruzioni in Zona Sismica 2007 



4 Le analisi di spinta 

Se si usa il metodo delle coordinate generalizzate, assumendo che il vettore spostamento sia 
descritto da una deformata <f), un vettore che da lo spostamento della struttura relativo a ciascun 
grado di liberta, modulato da uno scalare y(t), V equazione di equilibrio di una struttura a piu gradi 
di liberta: 

MX + KX = -MK g (1 
Pud essere riscritto come: 

M<S>y + K<S>y = -M\x g (2 

Che premoltiplicando per la trasposta della deformazione da luogo alia equazione (principio dei 
lavori virtuali): 

O r MOy + O r KOy = -O r MK g (3 

Si pud riscrivere la (3 nella forma consueta: 

my + ky = Lx (4 



Si supponga di aver normalizzato ad 1 la deformata in corrispondenza di un grado di liberta di 
piano. Con tale assunzione il valore di y coincide con lo spostamento del grado di liberta ove il 
vettore di forma vale 1. 

Dividendo per la massa generalizzata la (4 assume la usuale forma: 
y + a 2 y = px g (5 

Dividendo per il coefficiente di partecipazione p si ottiene: 

y 2 y (6 
—+»— = * 

p p 



La (6 pud essere vista come 1' equazione di equilibrio di un oscillatore semplice la cui risposta e 
data in termini di spostamento relativo: u=y/p. 
il + co 2 u = x\ (6 
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II termine L che nella (4 moltiplica I'accelerazione del terreno pud ess ere visto come la massa m 
di un oscillatore semplice equivalente: 

ft MI = m =L 



Moltiplicando primo e secondo membro dell (6 per m* di ottiene: 

(6' 



* y * 2 y 

m — Ym co — = m x 
P P 



Poiche: 
co 2 =k/m. 

sostituendo nella (6' si ottiene: 

* V * k V * .. 

m — Ym = m x 

p m p 



m — + (kp)— = m'x 



(6" 



(7 



II taglio alia base della struttura soggetta alle forze F=Kx =K(jjy, e dato da: V= F T 1= I T F. Se si 

usano le forme modali vale la relazione: 

F = K(/jy = ofMtpy (8 

Se si premoltiplica la (8 per il vettore I T di trascinamento orizzontale si ottiene, come appena visto, 
il taglio alia base della struttura: 

V{y) = I T F = I T K(/)y = I T co 2 M0y = affMfy = co 2 Ly = k—-y = k-p-y ^ 

m 

Pertanto sostituendo il risultato della (9 nella (7, si ottiene: 

..V(y) (10 
my + = Lx 

P 

Sino a che il legame e linear e si ha che V e proporzionale a y, pertanto V(y)/p=V(y/p): 
P P P P ■ 
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P P 

La (12 pud scriversi: 

y 

L-U + V{u) = Lx ;u = — 

P (13 

Se il legame e lineare tra V edy vale la (IS. Essa rappresenta la equazione di equilibrio dinamico 
dell 'oscillatore equivalente alia struttura. 

Quando I 'oscillatore semplice si sposta di u, la struttura si sposta diy= p u. 

La forza di richiamo ha il legame tra forza e spostamento della struttura. Pertanto la rigidezza k e 
proprio quella della struttura. Tuttavia, la forza massima alia base dell 'oscillatore equivalente non 
lineare si ha quando I 'oscillatore raggiunge uno spostamento u corrispondente a quello y di 
snervamento della struttura diviso per p. Ne segue che la massima reazione V che si pud avere 
nell 'oscillatore equivalente vale (V y )/p pari al taglio massimo di snervamento alia base della 
struttura diviso per p. 

Se si normalizza all 'autovettore in corrispondenza della copertura, la curva di spinta corrisponde 
a quella in cui in ascissa e riportato lo spostamento di sommitd (ove V autovettore vale 1) ed in 
ordinate il taglio alia base. 

La massa e data dalla somma delle masse pesate ciascuna per il valor e dell 'autovettore al livello 
della massa. 

II coefficiente di partecipazione p vale p= L/(<f)'M<f)). 
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4. 1 Caratteristiche dell'Oscillatore Equivalente ad una Struttura 



V 




Figura 11- Legame elasto plastico perfetto e curva reale 



Si consideri una struttura, che sottoposta ad un sistema di forze F distribuito in modo 
proporzionale al prodotto delle masse per un vettore di forma, abbia un legame taglio alia base 
spostamento in sommitd elastico perfettamente plastico, quale quello bilatero di figura. II il punto 
y y rappresenta lo spostamento in sommitd (se I'autovettore vale I in sommitd, sennd e il punto dove 
I'autovettore vale 1) cui corrisponde lo snervamento. 

II valor e d y * corrisponde a y y /p, il valore F y * corrisponde a V y /p. La massa deU'oscillatore e la 
somma delle masse pesate dall' ampiezza del modo in corrispondenza a ciascuna massa: 
(f) 1 MI = m* = L = y m ; • </> i 
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II coefficiente di partecipazione p e dato naturalmente da: 
0'MI _ m* 
jMi = P = J j m i tf 



F> 




Figura 12 legame delV oscillatore equivalente, ripreso dalla Circolare esplicativa delle Norme 
Italiane. 

Le curve di spinta reali sono difficilmente elastiche perfettamente plastiche. In genere hanno un 
andamento quale quello in Figura 12. Vi sono vari criteri per ottenere dalla curva non lineare 
reale quella equivalente elastica perfettamente plastica. 

La norma Italiana, nella circolare esplicativa, da il criterio riportato in Figura 12. Detta Fbu la 
resistenza massima del sistema strutturale reale ed Fb u *= FbJp la resistenza massima del sistema 
equivalente, il tratto elastico si individua imponendone il passaggio per il punto in cui il taglio alia 
base vale 0,6 F/, u , il punto di plasticizzazione (d y *; F y ), e quindi il valore F y si individua imponendo 
I 'uguaglianza delle aree sottese dalla curva bilinear e e dalla curva di capacitd per lo spostamento 
massimo d u * corrispondente ad una riduzione di resistenza <0,15 Fb u *. 

Si tratta naturalmente di un procedimento convenzionale, che pretenderebbe di dar luogo ad una 
struttura che da la stessa risposta in spostamento della struttura vera. In linea di massima la 
rigidezza iniziale: 

k*= F y / d y * 

e poco differente da quella elastica iniziale della struttura, ma comunque e minore. II periodo cosi 
determinato vale: 
T =2it(m /k ) 
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4.2 Valutazione delta risposta massima in spostamento dell'oscillatore 
equivalente. 

Vi possono essere vari modi per valutare la risposta dell ' oscillator e equivalente. Se si ha a 

disposizione un accelerogramma si pud valutare la risposta per integrazione al passo mediante 

analisi non lineare dinamica. L 'ipotesi di comportamento della struttura e elasto plastica. 

In realtd le strutture difficilmente hanno un comportamente elastico perfettamente plastico, 

pertanto la risposta va valutata con il "reale" modello non lineare ciclico della struttura. Esso 

potrebbe, in linea di principio, essere ricavato con una analisi di spinta ciclica. 

In generale si fa ricorso a modelli semplificati noti, quali appunto quello elasto plastico perfetto, 

oppure incrudente (se dopo lo snervamento la retta relativa al legame ha inclinazione diversa da 0, 

oppure un legame del tipo Tackeda, per strutture in cemento armato, o legami piu complessi. 

Un metodo semplificato per la stima della risposta non lineare, e quello di assumere che, se la 

struttura e a comportamento perfettamente plastico, lo spostamento massimo sia lo stesso della 

struttura elastica. 

. Figura 13 Criterio di uguaglianza tra massima energia di deformazione 



F 



u 




E el =l/2Kx u 2 

E p i=l/2Kx U p-l/2K (x up - x u ) 2 



y 




x, 



up 
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Questo criterio e valido per le strutture a periodo lungo. Per le strutture rigide in genere e piii 
approssimato il criterio di utilizzare Vuguaglianza dell'energia massima di deformazione tra 
struttura elastica ed elasto plastica. 
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Figura 14 legame tra spostamento massimo delta struttura elasto plastica ed elastica col criterio 
di uguale energia 



Si tratta di capire quando una struttura e flessibile (a periodo lungo) o rigida (a periodo corto). 
Fishinger [1999] ha suggerito di legare Vindividuazione delle strutture a periodo corto o lungo dal 
confronto con Vordinata T c dello spettro di risposta. E Vordinata ove finisce la zona di massima 
amplificazione della accelerazione (tratto costante dello spettro di progetto in accelerazione). 
Se T>Tc la risposta elasto plastica e quindi uguale a quella elastica, pari cioe alio spettro elastico 
in spostamento. 

Se T<Tc la risposta e data dalla seguente espressione: 



d = [l-(q-l)xi]-^l- 

L 1 Ji d(T*)xm* 
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Se il legame non e elastico perfettamente plastico, I ' espressione sopra indicata perde di validitd 
sempre piii al ridursi della dissipazione strutturale. Normalmente le strutture dissipoano meno 
rispetto al legame alasto plastico. Un criterio per capire la differenza e quella di sovrapporre il 
ciclo di isteresi della struttura a quello elasto plastico perfetto. II rapporto tra le energie dissipate, 
cie tra le areee incluse nei cicli da una misura della attendibilitd delle espressioni sopra riportate. 
Per ovviare a questo problema si possono utilizzare vari criteri. In alcuni casi sono state fatte 
analisi di confronto tra strutture che dissipano meno rispetto a quella elastica, trovando dei 
coefficienti di correzione (incremento) della risposta non linear e, che sono funzione del rapporto 
delle energie dissipate. 

Un criterio piu fisico, e quello di trovare la struttura elasto viscosa equivalente a quella elasto 
plastica. Si parte cioe dall'idea che se due strutture dissipano la stessa quantitd di energia debbono 
avere risposte dinamiche simili. Naturalmente questa ipotesi deve essere, ed in effetti, almeno sino 
a valori non eccessivi della energia dissipata, essa da risultati ragionevoli. Di questo problema 
tratta un notevole filone di ricerca noto appunto con il nome di "linearizzazione equivalente ". E ' 
possibile cercare la struttura lineare equivalente ottimale per ottenere diverse grandezze di 
risposta. Nel caso qui esaminato si cerca quella che da la stessa risposta massima in spostamento. 
Molti autori cercano invece la struttura che minimizza lo scarto quadratico medio tra la storia di 
spostamento della struttura non lineare e quella della struttura elasto viscosa equivalente. Se ad 
esempio sipesano maggiormente gli spostamenti massimi, i due obbiettivi tendono a coincidere. 



4.2. 1 Criterio di equivalenza tra struttura non lineare e struttura elasto viscosa. 

La storia della risposta di unsistema soggetto ad una forza: 

p(t) = p 0 sin (cot) E 1 

E del tipo: 

u(t) = u Q sin((i)t — 0) Eq. 2 

v(i) = u Q (i)sin((i)t — 0) Eq. 3 



Tralasciando la fase (j) I 'andamento temporale dello spostamento e della velocitd di risposta sono 
dati nella Figura 15. 
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038 0.75 1.13 15 



Figura 15 esempio andamento dello spostamento e delta velocita net tempo per eccitazione 
sinusoidale 

La forza dissipativa viscosa e data da: f(t) = cu(t), la forza elastica da: ku(t), la forza totale 
dalla somma delle due forze: /l(t) = cu(t) + ku(t). I diagrammi forza spostamento e forza totale 
spostamento sono indicati nella Figura 16 e nella Figura 17. 



f(t,T) 0 



-4189. - 5 




.4.641, 



fKUT) 



ku(t.T) 



-4.641.-5 




u(t,T) 



Figura 16 andamento delta forza viscosa con 
to spostamento 



Figura 17 andamento delta forza totale con to 
spostamento 



Se si fa un esperimento si misura in genere una parte elastica che si somma alia parte viscosa gid 
vista, ottenendo la curva di destra (Figura 17). L'energia dissipata vale 



2tt I a 



2n I co 



2k I co 



E D = j(cu)udt= jcu 2 dt= J c\om 0 cos(a>t — tp)Y dt = cncm 




poiche: c=2va> n m; m=k/(a>n) , lEq. 2'Eq. 4 diviene: 



E D = 2vco n k/(co n ) 2 ncouo = 2nv(co/co n )nkuo = [4nv(co/a) n )] f(l/2)k uq 2 ]; 



Eq. 5 
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E d = 47iv(co/co n )J E ei ; ove E d = [(l/2)k u 0 ]; 



Eq. 6 



Dalla Eg. 6 si ricava: 



v=[ 1/4 n] ( (oj a>) E d /E e! ; 



Eq. 7 



Dato un generico legame forza spostamento non lineare, si pud, per un ciclo di assegnata ampiezza 
ricavare i valori di rigidezza elastica secante K e dello smorzamento equivalente v. La prima si 
ottiene unendo Vorigine con il punto della curva ciclica forza spostamento corrispondente alio 
spostamento massimo, mentre lo smorzamento equivalente si ottiene con la Eq. 7. 




Figura 18. Esempio di legame el zsto plastico 



Se e data una curva elasto plastica incrudente, quale quella in 



Figura 18, se Q e V inter cetta della curva con I'asse delle ordinate, Ki la rigidezza elastica iniziale, 

K 2 la rigidezza del ramo plastico, d lo spostamento massimo, d y lo spostamento di snervamento, K 

la rigidezza secante, I 'energia dissipata pud essere espressa nella forma 

E d =4Q(d-d y ) ' Eq.8 
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Pertanto lo smorzamento si ottiene dalla espressione: 

v=[l/4n] (W co) E d /E e ,= v=[l/4n] (co n / co) 4Q(d- d y )/ f(l/2)k d 2 ] 



Eq. 9 



In genere si assume che il maggior contributo alio smorzamento complessivo e dato dalle frequenze 
prossime alia frequenza propria della struttura, quando (coj co)=l, pertanto la Eq. 9 diviene: 



v=[l/4n] 4Q(d- d y )/ f(l/2)k d 2 ]; 



Eq. 10 



Nel caso elasto plastico Q=F y , cioe I'intercetta e uguale alia forza di snervamento della struttura, 
nonche alia resistenza massima, e la rigidezza secante vale: k= F y /d, la Eq. 10 diviene: 



v=(2/n) (1- dy/d)= (2/n) (1- 1/ju); 



Eq. 11 




Figura 19 Smorzamento equivalente al crescere della duttilita in spostamento nelle strutture 
elasto plastiche perfette 

Si noti che lo smorzamento equivalente si addiziona alio smorzamento che comunque viene sempre 
associato alia struttura in campo elastico, in genere assunto pari al 5%. 

Quando la struttura ha un incrudimento positivo dopo lo snervamento, come indicato in 



Figura 18, 
kd=Q+k2d 

da questa si ottiene I 'espressione della rigidezza secante: 
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k=Q/d +k 2 . Eq. 13 
Inoltre, sempre con riferimento alia 

Figura 18, ed alia Eq. 10 si ha: 

d y =Q/(ki-k2); Q= kj dy-kidy Eq. 14 

detto: a= Q/ (dyki) Eq. 15 

v=(2/n)[a(fi-l)]/[( ju+a) juj Eq. 16 

Quando la derivata rispetto alia duttilita ju della Eq. 16 e nulla, si ottiene il valor e di ju per il quale 
lo smorzamento e massimo: 
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0.5 



jU ( vmax)=l+(l+a) 
vmax= (2a/ n) (l/[2(l+a) 1/2 +(2+a)J 



Eq. 17 
Eq. 18 



Si vede pertanto che, mentre nel caso elasto plastico perfetto lo smorzamento aumenta sempre 
all'aumentare della duttilita richiesta, nel caso di una struttura dotata di incrudimento lo 
smorzamento aumenta sino ad un certo valore della duttilita, oltre la quale lo smorzamento 
diminuisce.In realtd per le usuali struuture civili, il valore dell 'incrudimento e in generale solo di 
qualche percento. Se si utilizza VEq. 17 s ivede, che il valore della duttilita cui corrisponde lo 
smorzamento massimo e molto grande, per valori di dell 'icrudimento a=k2/kl <0.02. 




0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

a 



Figura 20 Duttilita a cui corrisponde lo smorzamento equivalente massimo al variare 
delV incrudimento 



Si pud pertanto concludere che sino a valori di incrudimento pari a 0.02 la dissipazione e quindi lo 
smorzamento equivalente aumentano con I 'ampiezza dell 'escursione plastica, cioe con la massima 
duttilita richiesta. Per valori di incrudimento maggiori, lo smorzamento cresce sino ad un certo 
valore della duttilita per poi diminuire. In particolare occorre controllare con attenzione i casi con 
incrudimento sinosuperiore a 0.1, valore per il quale la duttilita corrispondente alio smorzamento 
massimo e pari a 4.16, con uno smorzamento massimo pari a 0.33. 
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Figura Hsmorzamento equivalente alia struttura elasto plastica con incrudimento pari at 10% 
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5 Rappresentazione degli accelerogrammi in serie di fourier: lo spettro 
di potenza 

Dato un accelerogramma, ad esempio Tolmezzo NS del 1976: 
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Accelerogramma del Friuli 



Questo accelerogramma registrato ad intervalli di 0. 01 secondi per un totale di 3 728 punti, ha una 
durata di 37, 28 secondi. Pud essere approssimato con unafunzione del tipo serie di Fourier: 

n n 

f(t) = A n C0S + B n Sin <O n t = A 0 + Yj A, COS + B n sin 

0 1 

essendo : 
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In: 



il primo termine della serie e nullo in una funzione a media nulla quale un accelerogramma. 

Per esprimere I 'accelerogramma in serie difourier, occorre determinare il valore delle costanti An 
e Bn. Esse vanno scelte in modo che minimizzino lo scarto tra la funzione da approssimare aft), 
I 'accelerogramma, e la funzione approssimante fft). Se il numero di termini che si utilizzano e pari 
al numero di punti della funzione, lo scarto e evidentemente nullo, si potrebbe infatti scrivere un 
sistema di n equazioni in n incognite, ma il beneficio che si ottiene con la nuova formulazione e 
nullo. 

Se il numero di termini che si considera e modesto, si pud avere un beneficio reale, riducendo al 
dimensione del problema da trattare. 

In realtd, oggi non e tanto interessante questo aspetto, quante la comprensione di alcune alcune 
caratteristiche del dell 'accelerogramma che ci consentono di capire la fisica del problema. 

Le costanti An e Bn si ottengono cercando i valori che minimizzano lo scarto quadratico medio tra 
la aft) e la funzione approssimante fft): 

1 D 

S = -\(a(t)-f(t)fdt 

o 

derivado S rispetto ad A n eB n si ottengono le seguenti espressioni: 




il pregio della formulazione e che si pud inizialmente scegliere un numero limitato di armoniche 
approssimanti, valutare graficamente la differenza con la funzione da approssimare, e 
successivamente aggiungere un certo numero di termini sino a che il risultato non sia 
soddisfacente. Lo scarto diminuisce all'aumentare del numero di termini della serie. 

L 'espressione di Fourier pud ssere riscritta nella forma: 
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f(t) = f j C n (C0S(D n t-<l> n ) 

0 

C n =jA~TJ n 
<j) n = artang-^ 

Poiche si dimostra che, se N=cc, detta W la potenza media dell'accelerogramma: 

W = ±-]a(t) 2 dt = \f 4 C n 2 
u o 1 i 

un criterio per decidere come arrestare il calcolo di nuove armoniche e che la differenza tra la 
valutazione esatta della potenza e quella approssimata considerando un numero limitato di 
armoniche, sia modesta: 
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5. 1 Esercizio 

Dato il terremoto del Friuli 1976 (Tolmezzo N-S) 

Cliccando su QBasic si avvia il Basic. Si esegue quindi il programma ACC. 
II file accorciato a 16.64 secondi si trova su FR76C.TXT. 

Si esegua il calcolo senza AO con sole 10 armoniche. La maggiore e 2 lTlO/D= 62.8/16.64=3. 774 

Sullo schermo appaiono I 'accelerogramma di input, gli spettri dei Cn, funzione del periodo o della 
frequenza, il massimo valore f della funzione approssimante ed il massimo valore dello scarto 
quadratico medio tra funzione approssimante e da approssimare ( cioe I 'accelerogramma) . 

Si pud quindi chiedere il calcolo della risposta all 'accelerogramma. Viene prima calcolata la 
funzione di trasferimento e quindi data la funzione di risposta. 
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Si riporta 1' accelerogramma registrato a Tolmezzo (Friuli) nel 1976 in 
direzione N-S 
3728 0.01 
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